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The paper develops a mathematical theory of thyroid-pituitary inter- 
action. It is assumed that the pituitary gland produces thyrotropin, which 
activates an enzyme of the thyroid gland. The rate of production of thy- 
roid hormone is considered to be proportional to the concentration of that 
enzyme. It is further assumed that in the absence of the thyroid hormone 
the rate of production of thyrotropin is constant, but, in general, it is a 
linear function of the concentration of the thyroid hormone. This picture 
leads to a system of non-linear differential equations, which present 
great difficulties. This system, however, may be conveniently ‘‘linear- 
ized,’’ by considering that the relations between different variables are 
linear, but that within different ranges of the variables the coefficients 
are different. Using this approximation, it is possible to show that the 
system admits periodic solutions of the nature of relaxation oscillations. 

Such oscillations are actually observed in some mental disorders, 
such as periodic catatonia. The study of the effects of different para- 
meters of the system suggests different possible approaches to clinical 
treatment. In the light of this theory, the experimental determination of 
the parameters of the system becomes desirable and important. 


The system which controls the amount of thyroid hormone in the 
body is a feed-back regulator which can be described approxi- 
mately by a set of differential equations. In a previous paper, 
Danziger and Elmergreen (1954) presented the following non-linear 
equations in which @ and 7 are concentrations of system hormones 
at any time, ¢, and the remaining quantities are positive real 
constants: 

d0/dt = k,ma/(1 + m7) — 60 = P(6,7) (1) 
dn/dt =c —kn0/(1 +0) -g7 = Q(0,7) (2) 
Be ar, TOS 

This system of equations successfully describes most of the 

observed functions and malfunctions of the thyroid-pituitary ap- 


1 


2 LEWIS DANZIGER AND GEORGE L. ELMERGREEN 


paratus, but fails to account for the existence of sustained oscilla- 
tions of the hormone levels believed to produce the mental dis- 
order known as periodic relapsing catatonia. To show that sus- 
tained oscillations are not possible in the system of (1) and (2), 
we apply the First Theorem of Bendixon as given by Minorsky 
(1947), which states that limit cycles cannot exist in systems of 
this form if 0P/d0 + 0Q/d7 does not change sign. This negative 
criterion is satisfied by (1) and (2), and we can conclude that the 
nonlinearities present are not of the type which can produce sus- 
tained oscillation in this system. For purposes of stability study, 
therefore, the Langmuir adsorption terms in (1) and (2) may be 
linearized. 

To account for the observed sustained oscillations, we postu- 
late here a system of three first-order differential equations which 
are similar to the linearized form of (1) and (2), but which includes 
a specific mechanism by which the pituitary stimulates the thy- 
roid. It will be shown that this linear, third-order system will 
yield the same information as the set (1) and (2) and will also per- 
mit sustained oscillations of the hormone levels under proper 
conditions. 

The following symbols will be employed: 

7 is the concentration of thyrotropin at time, ¢; 

E is the concentration of activated enzyme at time, ¢; 

6 is the concentration of thyroid hormone at time, 4 

b, g and & are loss constants; 

a, h and & are constants expressing the sensitivity of the glands 

to stimulation or inhibition; 

c is the rate of production of thyrotropin in the absence of thy- 

roid inhibition. 

As in our previous communication, we assume that the rate of 
thyrotropin production is reduced by an amount proportional to the 
concentration of thyroid hormone and that the rate of loss of thy- 
rotropin is proportional to the thyrotropin concentration. Equating, 
then, the rate of change of 7 to the rate of production minus the 
rate of loss, yields 


dn/dt =(c-h0)—-ga7 when OS c/h (3) 
which is the linearized form of (2). The restriction placed on (3) 


is now required because the periodic solutions to be investigated 
will permit thyroid hormone levels of sufficient magnitude to sup- 
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press the production of thyrotropin. The production-rate term, 
(c — AO), can apply only for 6 < c/h as the pituitary gland can pro- 
duce no output in the presence of thyroid concentrations greater 
than c/h. To describe the system for large values of 0, we must 
include the degenerate form of (3), 


dn/dt=-gz when @2>c/h. (4) 


To describe the mechanism in the thyroid gland, we assume that 
thyrotropin activates a thyroid enzyme, which, when activated, 
catalyzes the production of hormone, 6. Such an effect of thyro- 
tropin has been observed by Vanderlaan and Greer (1950), by Ghosh, 
Woodbury and Sayers (1951), and by others, who reported that the 
thyroid of hypophysectomized animals was unable to trap radioac- 
tive iodide, and that this ability was restored by the administra- 
tion of thyrotropin. While it is not necessary to identify the en- 
zyme involved, we may imagine it to be a peroxidase which oxi- 
dizes iodide to iodine before its incorporation into the thyroid 
hormone molecule. Thyroid hormone production, according to this 
assumption, will depend on the concentration of activated enzyme 
and not directly on the level of thyrotropin. To describe this en- 
zyme action, we equate the rate of change of activated enzyme to 
the rate of activation, assumed proportional to 7, minus the rate of 
inactivation, assumed proportional to E. This yields 


dE/dt = mr —kE. (5) 


Similarly, equating the rate of change of thyroid hormone concen- 
tration to the rate of production minus the rate of loss, yields 


d0/dt = aE — b6, (6) 


which is similar to the linearized form of (1) except that the 
production-rate term is now a function of the activated enzyme 
level rather than the thyrotropin concentration. 

Equations (3), (4), (5), and (6) now represent the feed-back regu- 
lator action of the thyroid-pituitary system. The physical condi- 
tion that the hormones and enzyme cannot have negative values is 
now inherent in these equations and need not be stated separately. 

This system agrees with the observations in various clinical 
situations. If the pituitary gland is removed, ¢ will become zero, 
7 will fall exponentially with time, £ will follow z, and @ will fol- 
low E. If E is unable to function, because of the presence of some 
enzyme poison or of inadequate supply of iodide, the rate of pro- 
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duction of thyroid hormone is reduced in magnitude, and the re- 
sulting increase in 7 will not overcome the effect of the poison or 
the deficiency. If the thyroid gland is removed, @ will fall ex- 
ponentially with time, and 7 will approach the value ¢/g. 

A more convenient form of the equations may be obtained by 
means of variable changes. Let 
a2 =(g/h)7 be a new variable proportional to 7; 
y =(gk/hm)E be a new variable proportional to E; 
C =c/h be a controlling parameter; 
K =(ahm/bgk) be a dimensionless overeall gain constant; 
T, = 1/9, T, = 1/k and T, = 1/6 be system time constants. 
Substitution of these quantities in (3), (4), (5), and (6) yields 


T,da/dt+a=C-@ when @SC (7) 
T,da/dt +#=O when @62C (8) 
T,dy/dt+y=2 (9) 
T,d0/dt +0 = Ky. (10) 


As a first step in the analysis of the set of equations (7) through 
(10), the steady-state levels of the dependent variables will be de- 
termined. For a normal system, the steady state may be charac- 
terized by zero rates of change of the hormone and enzyme levels; 
further, the steady-state level of 6 must be less than C. Setting, 
then, the derivatives in (7), (9), and (10) equal to zero and solving 
for the steady-state levels of the hormones yields 


x, = C/(1 + K) (11) 
0, = CK/(1 +k). (12) 


These steady-state equations correlate the behavior of the thyroid- 
pituitary regulator with the system constants. The parameter 
C =c/h is seen to control the steady-state levels of both hor- 
mones, while the gain constant K =(ahm/bgk) has a large effect 
on #, and a relatively small effect on the steady-state level of thy- 
roid hormone. It can be shown that typical values of K will be 
sufficiently large to permit the approximation K/(1 + K) = 1 so that 
(12) may be written as 


0% C =o/h. (13) 


Since the controlled quantity in the thyroid-pituitary system is 
the level of thyroid hormone, equations (12) and (13) are of primary 
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importance. These equations show that if the thyroid gland-is 
normal the steady-state level of thyroid hormone is essentially de- 
pendent on the pituitary constants c and h, whereas if K is large 
the effects of the sensitivity constants a and m and of the loss 
constants 6, g, and & are almost completely nullified by feed-back 
action. This result is analogous to that of feed-back amplifiers or 
servo-mechanism systems where the output-input ratio can be made 
practically independent. of system parameters. The thyroid- 
pituitary system, therefore, can be considered to be a feed-back 
regulator which, in terms of the servo-mechanism analog, has an 
input C, an output 6 and an open loop gain K. 

System malfunction characterized by abnormal steady-state hor- 
mone levels, such as hyperthyroidism or hypothyroidism, in con- 
ditions not involving deficiencies of raw materials or enzymes, 
arise, as shown by (13), from a defect in the pituitary gland which 
would produce an abnormal ratio c/A. In toxic goiter, for example, 
the presence of excessive levels of thyrotropin, as well as of thy- 
roid hormone, indicates a higher than normal c/h. Thyroidectomy 
would reduce the level of thyroid hormone but would increase the 
level of thyrotropin. This can be seen by making the constant, a, 
in (6) and hence the value of K, zero. From (11) and (12), then, 
a becomes equal to C and @, becomes zero. From this we con- 
clude, in accordance with clinical observation, that thyroidectomy 
would relieve some, but not necessarily all, of the symptoms of 
toxic goiter. Thyroidectomy followed by administration of suf- 
ficient thyroid extract to keep the thyrotropin concentration low 
would seem to be a more efficient treatment of the disorder. 

Hypothyroidism may result from a variety of causes. The pitui- 
tary type is due to abnormally low values of c/h and, hence, to a 
malfunction of the production mechanism of the pituitary gland; it 
is associated with low values of thyrotropin. The type caused by 
goitrogens, which inhibit EZ, is due to failure of the production 
mechanism of the thyroid gland; it is associated with high values 
of thyrotropin. Analysis of the steady-state levels for the condi- 
tion of iodine deficiency requires further discussion of the produc- 
tion-rate term in (6), which has been assumed to be linear for all 
positive values of activated enzyme, E. This term will be linear 
only as long as a sufficient amount of iodine is present to permit 
thyroid-hormone production at the rate demanded by the level of EF, 
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If iodine intake is severly limited, the term ak will express the 
thyroid production rate only so long as the rate of iodine utiliza- 
tion does not exceed the rate of iodine intake. Mathematically, 
we may define a critical value of enzyme, EF, for which iodine is 
used as fast as it becomes available. Then for FE > E_, the thyroid 
production rate would be constant at the value af. If the normal 
steady-state level of activated enzyme is greater than E, the sys- 
tem is limited by the supply of iodine, and the steady-state thyroid 
hormone level must be less than normal. 

The last system malfunction which we shall discuss is the men- 
tal disorder called periodic relapsing catatonia (Gjessing, 1939, 
1953), which is characterized by fairly regular variations in the 
basal metabolic rate (BMR) and in the severity of the psychosis. 
In applying equations (7) through (10) to this condition, we assume 
the observed variations in the BMR are due to sustained oscilla- 
tions in the thyroid hormone level and that the periodic nature of 
the solution of the system equations is due to instability for 6<C, 

To investigate the stability of the system, we obtain the third- 
order differential equations describing the variation of thyroid hor- 
mone with time. Eliminating x and y from (7), (8), (9), and (10) 
yields the pair of linear equations 


a,d°0/dt* + a,d°0/dt® + a,d0/dt+(1+K)@=KC, @<C (14) 


a,d°0/dt® + a,d°0/dt? +a,d0/dt+0=0, O>C (15) 
where the coefficients are the positive real numbers given by 

a, =T,T,T, 

agel Tat +27, (16) 


a,=T,+T,+T, 


Employing the classical operator D = d/dt, the characteristic 
equations of both (14) and (15) may be written as 


a,D°+a,D*+a,D+a, =0, (17) 


where a, = 1+K for (14) and a, = 1 for (15). 
Stability of (14) and (15) may be discussed in terms of a face 
tored form of (17) which is 


[D — BIL(D - a)? + w?] ‘ (18) 
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where &, 6 and o are functions of the coefficients. The roots of 
(17), expressed in terms of the factors of (18) are 


Depo p= ta, (19) 
where @ may be real or imaginary. 


A general solution of both (14) and (15) may now be written in 
the form 


0=0 +Ne*sM ce sin (wt + ), (20) 


where N, M and © are determined by the initial conditions. Clearly, 
all six of these constants differ for the two solutions included in 
(20): for (14), 0, = CK/(1+ K) and the solution is valid for 6 < C; 
while for (15), 6 =0 and the solution is valid for 6>C. Further, 
a, differs for the two. cases and hence a, 8 and @ will necessarily 
-have different values. Stability of (14) and (15) requires that 


lim @= 6 which obtains only if both « and 8 are negative. Anal- 
t—?> co 


ysis of (17) and (18) shows that B must always be negative for 
positive real coefficients, while « may take on either positive or 
negative values depending upon the relative values of the coeffi- 
cients. A stability criterion for (14) and (15) may then be developed 
by finding a relation connecting the magnitudes of the coefficients 
which will encompass all possible values for which «<0. This 
can be done by considering the critical case with « =0 which 
separates the stable solutions with «<0 from the unstable with 
a>0. The quadradic factor in (18) becomes (D? + w”), for a= 0, 
and can be divided out of (17) leaving a linear factor (a,D + a,) 


and a remainder (a, — a,@”) D +(a - ao”). If (D? + w”) is to be 
a factor of (17), both terms in the remainder must vanish, yielding 
wo? = a,/d, = a/a, (21) 


and 


a =a, Ay/ Gy ; (22) 


The critical solution is then defined by a=0, B =~a,/a, and 
@* = a,/a,, and occurs when the coefficients are related by (22). 


In order to employ (22) as a stability criterion, we need only 
observe that if a =0, the characteristic equation becomes one of 
second degree which for positive real coefficients must always 
have roots with negative real parts. Therefore for small a, the 
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third-order differential equations (14) and (15) are stable with 
a <0; while for a, large enough, a > 0 and instability will exist. 
This leads logically to the criterion that (14) and (15) are stable 
for all values of coefficients such that a, <a,a@,/a,. This result 
can also be obtained by employing any of the stability tests appli- 
cable to linear differential equations with constant coefficients. 
See, for example, the criterion of Routh or Hurwitz as given by 
MacMillan (1951). 

To apply this stability criterion to (14) and (15), we employ (16) 
to express the coefficients in terms of the time constants. As 
T,, T, and T, are positive real numbers, it can be shown that 
@,a,/a, has a minimum value of nine when the time constants are 
equal. Hence, for (14) in which a, = 1+ K, the solution (20) will 
be stable for any value of K less than a critical value defined as 


K~ = 4,0,/a, —<1, (23) 


where K is greater than or equal to 8 and corresponds to the criti- 
cal solution with a=0. For (15), in which a, =1, stability will 
exist for all values of the time constants, the roots of its charac- 
teristic equation being —-1/T,, ~1/T, and -—1/T, and, hence, all 
real and negative. 

Consider now the case where K >K,. For any set of initial con- 
ditions such that 9 <C, an unstable system governed by (14) will 
exist with a solution of the form of (20) in which a>0. This will 
yield an oscillatory variation of @ with the amplitude of the peri- 
odic term increasing exponentially with time. Ultimately, in this 
mode of variation, @ will reach the value C, and the system will 
degenerate to the stable system of (15) with a zero steady-state 
level. To show that the stable solution must exist for a finite time 
interval, we may use (10) which shows that d0@/dt must be con- 
tinuous as neither @ nor y may change value instantaneously. Thus, 
as d6@/dt is positive when 6=C in the unstable mode, the thyroid 
level must continue to rise before it can begin to fall toward the 
zero steady-state level of (15). This aperiodic variation of the 
thyroid hormone can continue only so long as @>C, and when @ 
falls to the value C, (14) will again govern the system. Here, 
again, the continuity of d6/dt requires that @ fall below the level 
C before rising in the unstable mode. Clearly, this sequence of 
events will give rise to a sustained oscillation of the thyroid hor- 
mone concentration about the value C, governed alternately by (14) 
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and (15) and would obtain for any set of initial conditions. This 
type of solution, caused by periodic degeneration of an unstable 
System, is a relaxation oscillation (Minorsky, loc. cit.). 

Any treatment of periodic relapsing catatonia must suppress the 
oscillation of the hormone levels, which has been shown to be 
caused by abnormally large gain, K =(ahm/bgk). As it is not fea- 
sible to change, by treatment, any of the system constants so as 
to reduce K, the only procedure which will abolish these oscilla- 
tions is one which would suppress the thyrotropin output of the 
pituitary gland and the output of the thyroid gland. ‘The level of 
thyroid hormone necessary to maintain the BMR would then have to 
be supplied from an external source. Thyroidectomy followed by 
administration of thyroid extract would meet these conditions if 
the thyroid extract added to the system were sufficient to limit the 
thyrotropin level to safe magnitudes. The administration of suf- 
ficient thyroid extract will, however, produce the same results 
without thyroidectomy. This procedure will be shown to comprise 
the proper treatment of periodic relapsing catatonia. 

If thyroid extract is administered in daily doses such that the 
effective rate of assimilation is a constant, FP, the original equa- 
tions must be modified to include this second source of thyroid 
hormone. For this modified system, equation (6) becomes 


d0/dt =ak +R -b0 (24) 
and (10) becomes 
T,d0/dt+0=Ky+T,h. (25) 


The total system thyroid, 6, can now be considered to consist of 
two components: the endogenous hormone, z, and the exogenous 
hormone 7,f. Substituting 6=2+T7,F in (7), (8), (9), and (25) 
with C’ = C - T,R, yields 


Tde/dt+*=C’-~2 when 25C’20 (26) 
T ,dx/dt + “=O when 22C’20 (27) 
T,dy/dt+y=« (28) 
T,dz/dt+2=Ky (29) 


for the system under treatment. Note that C’ is proportional to the 
rate of thyrotropin production in the absence of endogenous thyroid 
hormone inhibition and cannot be negative. If & 2 C/T,, C’ must 
be zero, and (26) becomes (27) for all values of z. The behavior 
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of this modified system, therefore, will differ for the two cases: (a) 
R <C/T, and (b) Rk > C/T;,. 

In case (a), where the rate of assimilation of exogenous thyroid 
hormone is small, (26) through (29) apply. These equations are the 
same as the original set (7) through (10) except that 2 replaces 
6 and C’ replaces C. Considering the endogenous hormone con- 
centrations, # and z, the administration of thyroid extract such that 
R <C/T, is equivalent to reducing the controlling parameter C to 
C’, and all of the preceding analyses and conclusions for the un- 
modified system are applicable. The essential difference between 
this modified system and the original is the effect on the total sys- 
tem thyroid hormone level. Equation (12) will apply for the steady- 
state endogenous thyroid hormone with C = C’ but the steady state 
of the total system hormone, being the sum of z, and T,k, would 
become 

6. =(KC +T,R)/(1+K), (30) 
which is slightly higher than that of the unmodified system. Thy- 
roid extract treatment, therefore, would not be effective in toxic 
goiter as the thyroid level would increase slightly, although some 
benefit might result from the reduction of the thyrotropin level. Im- 
provement would result in cases of hypothyroidism with thyroid ex- 
tract treatment regardless of the cause of thyroid deficiency. In 
these conditions administration of thyroid extract sufficient to 
bring the thyroid hormone level up to the required steady-state 
value would reduce the steady-state thyrotropin concentration with- 
out harmful effect. In periodic relapsing catatonia, if R < C/T, , no 
improvement will result. In this case, oscillations of 2 about C’ 
will exist, yielding, as in the original system, oscillations of 0 
about C; the only difference will be a reduction in the average 
values of the periodic thyrotropin and enzyme variations. 

In case (b) where the rate of assimilation of exogenous thyroid 
hormone is sufficiently large to suppress the production of thy- 
rotropin, only equations (27), (28), and (29) apply. Here R > C/T,, 
and the system, for any K and any set of initial conditions, is 
absolutely stable. For unequal time constants, the following gen- 
eral solutions apply: 

a= X, en t/Ts (31) 
pel ei eat Te (32) 


ga Z, eV 4 Zee ee Leer (33) 
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where ¢=0 at the instant when treatment starts; X, = %, Y, and 
Y, are determined by @ and y%; and Z,, Z,, and Z, are determined 
by 2%), Yo, and 2,. These aperiodic transient solutions have zero 
steady-state values and are independent of the magnitude of R; 
they are completely determined by the time constants, the gain, K, 
and the state of the system at t=0. The system with R > C/T, 
would ultimately become quiescent with the thyroid hormone being 
supplied entirely by the external source. For treatment of periodic 
relapsing catatonia, the rate of intake of exogenous thyroid hor- 
mone should be & = C/T, so that a steady state, 6.=C, is 
achieved. This results in complete stabilization of the thyroid 
hormone concentration at a value which is only (K + 1)/K greater 
than the normal steady state of (12). If C is abnormally small, then 
FR should be sufficiently large to raise 6. to the level required for 
health. 

The transient variation of endogenous thyroid hormone, 2, will 
determine the overshoot of the thyroid hormone concentration fol- 
lowing the start of treatment with large F,. As 0=2+T,F and as 
the steady state is 0. = T,h, the single maximum of 2, from (33), is 
the overshoot. This peak value of hormone can reach dangerous 
magnitudes, and care must be taken in the administration of thyroid 
extract to hold the overshoot to a low value. Several possible 
treatment procedures which will minimize the overshoot in thyroid 
hormone are known. Gjessing (Joc. cit.) originally favored starting 
treatment with large values of F and suggested, as an optimum 
starting point, a time when the BMR is near a minimum in the pre- 
treatment oscillation. This starting point was determined by ex- 
periment and produced satisfactory results. We assume that vari- 
ations in the thyroid hormone level are reflected in variations in 
the BMR; and, hence, Gjessing’s method is equivalent to starting 
treatment when the oscillatory thyroid hormone is a minimum, ‘Since 
- the thyroid hormone level will jump by the amount 7,F soon after 
treatment is started, it is evident that a minimum initial value will 
obtain if 6 is a minimum at t= 0. The endogenous hormone, 2, from 
(33) will have a maximum after ¢ = 0, however, and a more precise 
identification of the optimum starting point is possible. This would 
depend upon the system parameters as well as the initial condi- 
tions but would not be of much use clinically as none of these 
quantities can be accurately determined. 
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A second procedure, suggested by Danziger and Kindwall (1953), 
starts with small FP at any point in the oscillation cycle. The daily 
intake of exogenous thyroid hormone is then gradually increased 
until a stabilizing level is achieved. This method permits initia- 
tion of. treatment without a long period of preliminary study and 
will yield safe peak values of thyroid hormone. 

A third possibility is suggested by inspection of (25), which 
shows that a maximum @ results when d@/dt = 0 and y has some 
value, y,- This maximum is then 


6. Ky, + Tk, (34) 


max 
where the overshoot is Ky,. If the thyroid enzyme be inhibited 
before treatment starts so that y, is small, treatment with large Rk 
could be safely started without regard to the state of the system. 
This result could be achieved by administration of a thyroid enzyme 
poison before the start of thyroid extract treatment. 


General Remarks 


The preceding analysis illustrates the utility of mathematical 
representation of an endocrine control system. The feed-back 
aspects of the thyroid-pituitary homeostatic mechanism make pos- 
sible an analogy with this and perhaps other physiological phenom- 
ena on the one hand, and feed-back amplifiers or servo-mechanism 
systems on the other. A major difference in the philosophy of study 
in these fields arises from the fact that the electronic or servo- 
mechanism engineer is interested primarily in synthesis whereas 
the endocrinologist must necessarily consider analysis of a system 
in which measurements are difficult and numerical values of para- 
meters are either unknown or known with relatively little accuracy. 

In the thyroid-pituitary homeostatic mechanism, particular solu- 
tions of the equations presented are of secondary importance. If 
the system parameters are known or can be estimated, standard 
procedures can be used to obtain particular solutions. We employed 
high-speed electronic analogue computer techniques to obtain rapid 
checks on the conclusions made. For this computer study system 
parameters were estimated, partly from the literature and partly 
from the theory developed. This work, Danziger and Elmergreen 
(1956), illustrated additional advantages of mathematical analysis 
in the determination of some of the parameter values. In this case, 
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the time constants, T, and T,, could be estimated from data in the 
literature while the time constant, T,, of the thyroid enzyme is 
unknown. Using relaxation oscillation theory of periodic relapsing 
catatonia and clinically observed periods of oscillation, the time 
constant 7, was evaluated in terms of 7, and 7,. Further, an esti- 
mated value of typical gain, K, was obtained by computing the 
critical value from (23). The theory developed provided a means 
for estimating, with a minimum of experimental data, the numerical 
values of parameters in a typical system. 


The expenses of this investigation were defrayed, in part, by 
special gifts of Miss Margaret Kckart, Oak Park, Illinois, and the 
Willard Tablet Company, Chicago, Illinois. 
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Given a two-person continuous non-constant sum game, a ‘‘solution’? 
can be arrived at even in ignorance of the pay-off functions as each 
player tries to maximize his own pay-off by trial and error. The solution 
of the game considered here is either the intersection of two ‘‘optimal 
lines’? (the case of the stable equilibrium) or parasitism of one player 
upon the other (the case of unstable equilibrium). The introduction of 
secondary satisfactions (linear combinations of the pay-offs) affects both 
the position and the stability of the equilibrium and thus the resulting 
primary satisfactions (the pay-offs). These effects of the matrix of trans- 
formation are investigated. In particular, it is shown that the social 
optimum, which is also the solution of the game where collusion is al- 
lowed, is attained if and only if the columns of the matrix are identical, 
and that this solution is generally stable. Some suggested connections 
with biological situations are discussed. 


In the work of N. Rashevsky (1951) and in the author’s (Rapoport 
1947a,b), situations are discussed in which each of two individu- 
als is faced with a decision to choose one of a continuum of val- 
ues designating an ‘‘effort.”” The pair of values thus chosen maps 
upon a pair of “‘satisfactions’’ through two specified functions per- 
taining to the two individuals respectively. The problems deal 
with the determination of the efforts as both individuals attempt to 
maximize their respective satisfactions. 

Usually the efforts are determined by the intersection of the two 
‘‘optimal lines,”’ 

dS ,/dx#=0; AS,/dy =0, (1) 


where x and y are the respective efforts, while S,(z, y) and S,(a, y) 
are their satisfaction functions. 

The situation is typical of the two-person game in normal form. 
As is well known for the case of the zero-sum game, where S, = 
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-S,, a solution is the so-called minimax solution in the space of 
mixed strategies, if it is defined as the best that each player can 
do for himself when the functions S, (there called the pay-off func- 
tions) are known to both. The theory of the minimax solution does 
not apply non-equivocally to the general two-person game because 
it rests on the assumption that the interests of the individuals 
(players) are opposed at all times—what one gains, the other loses. 

In subsequent developments of game theory, a solution analogous 
to that determined by equations (1) has been called a Nash equi- 
librium (Nash, 1951). Such an equilibrium is characterized by the 
fact that neither player wishes to depart from it by his own effort, 
which is the only variable he controls, since departure in either 
direction results in immediate decrease of his own satisfaction. 

The Nash equilibrium, however, as it is developed in game the- 
ory is not determined by the intersection of two ‘‘optimal lines,”’ 
given by (1). There the continuous game derived from a finite game 
by the introduction of mixed strategies results in pay-off functions 
which are linear in the variables representing the choice of mixed 
strategy. Hence the maxima of those manifolds cannot be obtained 
by setting derivatives equal to zero. Rather, the maxima appear as 
a result of confining the variables to the intervals (0, 1) since they 
are essentially probabilities. Thus the theory revolves, as it does 
to some extent in linear programming, around considerations of the 
properties of convex spaces and other matters rather removed from 
classical equilibrium theories. 

In this paper, we shall be concerned with the classical equilib- 
rium, which is a Nash equilibrium only if one takes as one’s model 
not a finite but a continuous game and forgoes the consideration of 
mixed strategies, which in the case of a continuous game involve 
function spaces, Such an equilibrium is also central in Rashevsky’s 
treatment of hedonistic behavior (op. cit.). We shall therefore refer 
to it as the Rashevsky-Nash equilibrium. 

The Rashevsky-Nash equilibrium is interesting because, in the 
process of the players’ seeking maximum satisfactions, it can be 
arrived at even if the players are ignorant of the functions S, and 
S,. This happens, however, only if the equilibrium is stable (see 
below). Whether it is stable or not, however, the equilibrium need 
not be the best that each player can do. If S, and S, are known and 
if the players can engage in collusion, they can find points (x, y) 
which result in greater S, and S, than the values determined by the 
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equilibrium. They can, for example, choose an (z, y) so as to max- 
imize §,+8,. But if they do so, problems arise about-how to 
divide the spoils, i.e., problems of bargaining. The problems 
become interesting and difficult if the satisfactions are non- 
transferable (being, for example, of subjective nature) or non- 
conservative. Other problems arise if the players for want of pos- 
sibility of communication or for other reasons cannot engage in 
collusion. In the latter case, a very interesting situation can arise 
where both players act “‘rationally’’ and yet the ‘‘social optimum” 
solution, where S, + S, is maximized, fails to obtain. Consider the 
game whose pay-offs are indicated in the following matrix (Raiffa, 
1951):* 


B, B, 
Ol, ek +9) (~10, +10) (2) 
Gel(e10N= 10) (0, ct 


Here player I is to choose a row, player II a column, and the num- 
bers in parentheses indicate their respective pay-offs. Player I 
rightly concludes that he will be better off in choosing %,, no mat- 
ter what player II does. Accordingly, in view of the ‘“‘sure thing 
principle’ (Savage, 1954), he chooses a. Player II argues ex- 
actly the same way, and as a result, the entry (a, 8,) = (-9, -9) 
is chosen where both players lose, whereas the ‘‘social optimum”’ 
is (+9, +9). 

In this paper, we shall be concerned only with the solutions de- 
termined by the Rashevsky-Nash equilibria because it seems that 
these solutions are more applicable to biological situations. Typi- 
cally biological adaptations in an organism (as distinguished from 
the species) are reflections of what does happen rather than of what 
may happen. A homeostatic device, for example, insures the per- 
petuation of a Jocal optimum through the control of, typically, a 
single variable rather than a search for a general optimum based on 
the knowledge of the over-all situation and the exchange of this 
knowledge. 

The interdependence of homeostatic devices is not precluded 
however. Such interdependence is equivalent to communication, 
but is typical of communication about ‘“‘here and now’’ and does not 


*The game is known also in earlier literature under the nickname, 
‘*The Prisoner’s Dilemma.”’ 
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take cognizance of the entire situation. It may, however, give an 
appearance of ‘‘teleological’’ adaptation and thus provides a valu- 
able theoretical basis for ‘‘mechanical’’ explanations of teleologi- 
cal effects. We shall introduce such interdependence between the 
satisfaction functions of the players to show how the equilibrium, 
the only attainable solution, is thereby modified both in position 
and stability. In special cases, the social optimum can be made to 
coincide with the equilibrium. 

We call S, and S, the primary satisfactions and introduce sec- 
ondary satisfactions, 


S =48,+B8, (3) 
Sam Cs, DS, 
where A, B, C, D, are independent parameters. The matrix 


AB 


M (4) 


C D 


will be called the social matrix. We shall examine the equilibria 
pertaining to S, with regard to the resulting primary satisfactions 
S;; hence, the effect upon S, of the matrix M. 

The matrix M can be given various interpretations. In the work 
of N. Rashevsky (op. cit.), it describes the ‘‘altruism’’ prevalent 
in a social system. Thus if M is diagonal, the society is ‘‘ego- 
tistic’’ since each member considers in his secondary satisfaction 
only his own primary satisfaction. If the entries of M are all equal, 
the society is ‘‘egalitarian,’’ all for one and one for all. If M has 
only one column of positive elements, the society is a despotism, 
oriented toward the satisfaction of a single member. If the rows 
are all monotone decreasing, the society is hierarchical, etc. Neg- 
ative entries represent ‘‘antipathies.”’ 

The effect of this matrix on the position of the equilibrium was 
discussed by Rashevsky (op. cit.). He did not, however, discuss 
the effect of M on the stability of the equilibrium, an important con- 
sideration because for the unstable case, as we shall see, the be- 
havior of the individuals is not represented by the equilibrium but 
degenerates, at least in the game here considered, into parasitism, 
where one of the players brings his effort to zero. Moreover, in 
discussing the change of position of the equilibrium, Rashevsky 
does not indicate the resulting change in primary satisfactions. 
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As we shall see, we may have, as a result of this effect of the 
matrix M, an increase in the efforts expended but a decrease in the 
associated primary satisfactions. 

Questions concerning M have been raised also in game theory. 
If M is diagonal but not scalar, it represents a set of independent 
homogeneous linear transformations on the set of pay-offs. In J. 
von Neumann’s theory of utility (von Neumann and Morgenstern, 
1947), the utility function of an individual is defined only up to a 
linear transformation. In this context where utility is defined not 
in terms of a universal unit like money, but in subjective terms, it 
is impossible to determine the unit and the origin of an individu- 
al’s utility scale. Therefore, it is desirable to have conclusions 
of game-theoretical theorems invariant with respect to independent 
linear transformations of the utility functions involved. Nash and 
Raiffa (op. cit.) have studied classes of games which exhibit such 
invariance. As it turns out, Nash equilibria are invariant with re- 
spect to independent linear transformations on the separate utili- 
ties, but not with respect to the general linear transformations on 
the vector of utilities. In other words, Nash equilibria are neces- 
sarily invariant only if M is diagonal. But we shall consider the 
effects of a more general class of M. 

In biological situations, M may have the following significance. 
Let two organisms produce substances in the amounts 2 and y re- 
spectively and exchange portions of them. Let the efficiency of 
some biological function depend on the amounts of the substances 
apportioned to each individual and also on the efforts expended in 
producing the substances. Then, if each organism through its ho- 
meostatic devices seeks to maximize his efficiency S,, we have an 
analogue of the two-person game. As will appear below, given cer- 
tain forms of S,, if the Rashevsky-Nash equilibrium is unstable, 
one organism will become ‘“‘parasitic’’ on the other, i.e., cease to 
produce his substance. If, however, the efficiencies are ‘*hooked 
up’? by means of the social matrix, the situation will, in general, 
change the position and sometimes the stability of the equilibrium. 
The question naturally arises whether in this case we should not 
speak of two interdependent homeostatic devices of a single or- 
ganism instead of two separate organisms. The question is a mat- 
ter of definition. We may define parasitism as the case where one 
organism does not contribute to the common stock; symbiosis as 
the case where both contribute, but the ‘‘satisfactions”’ are not 


20 ANATOL RAPOPORT 


interlocked, i.e., neither responds to the needs of the other; inte- 
gration as the case where the satisfactions are interlocked. One 
can thus speak of varying degrees of integration, and indeed, in 
biology the line between symbiosis and integration is not sharply 
drawn.—It is just possible that the evolution of the integrated sin- 
gle cell organism from symbiotically related simpler virus-like, or 
gene-like organisms was a process of progressive integration. 

We shall confine our attention to a particular two-person non- 
constant-sum game in normal form, which is a slight generalization 
of the situation treated in a previous paper (Rapoport, 1947a, here- 
after called I). The game is defined by the equations 

S, = log (1 + px + gy) - Bz (5) 

S, = log (1 + ga + py) - By, (6) 
where 20, y20, p, g, and B are positive constants, p+ q=1; 
B <p. The last inequality is imposed to avoid the consideration 
of negative §,. 

In other words, each player exchanges a fixed fraction gq of his 
own product for an equal fraction of the other’s. The satisfactions 
depend positively on the amounts received and negatively on the 
efforts, where B defines the unit of effort. The optimal lines are 
obtained by equation 


P 


OS / Oa = nen -nne teen 26 7 
1+ pe+gy (7) 

0S,/dy = a----2--- Le ees -B =0. (8) 
1+ gu + py 


If p 4 q, they are two straight lines, and the equilibrium is their 
intersection 


G=7=p/B-1. (9) 

An equilibrium will be called stable if, in its neighborhood, the 
point determined by the efforts of the players tends to move toward 
the equilibrium as the players attempt to maximize their own re- 
spective satisfactions. If the point moves away from the equilib- 
rium, it is called unstable. 

Theorem 1. Let y,(x) and y,(x) be the equations of any two opti- 
mal lines. Then the equilibrium is stable if and only if in its 
neighborhood |y; 

Note: We conceive of the motion of the point (2, y) as follows. 
As one coordinate remains stationary, the player who controls the 


| 
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other brings the point (2, y) into his own optimal line. Then his 
coordinate remains stationary, and the other player does the same. 
If we think of the time intervals during which one coordinate re- 
mains stationary approaching zero, the point (z, ¥) will describe a 
smooth path in the z—y plane. It is easy to see that in the neigh- 
borhood of intersection of the optimal lines, this path will move 
either toward the intersection or away from it. 

Proof of the theorem. Let player I choose some value, 2‘!, near 
(0, 0). Player II will choose a value y™ such that the point(«™, 
y™) is on the curve y, (x); in other words, the first pair of choices 
result in [2™, y, (2)]. Now player I again chooses a value a” so 
as to bring the point to his own optimal line, y, (v). This value is 
the intersection of the horizontal line y = y,(«™) with y,(z), which 
is approximately 


(@)) ° 


y;(0) —-y, (0) 


We shall then have |2™| < |x | or |x” | >|2™| depending on whether 
y(0)| <ly{(0)| or |y%(0)|>]y{(0)|. Hence the above inequalities 
determine whether in the neighborhood of an equilibrium player I’s 
successive choices approach the equilibrium or depart from it, and 
the theorem is proved. 

Let us now examine the equilibrium determined by the equations 
(7) and (8). Rewriting the equations as 


pe+qy=p/B-1, (11) 
qz+py=p/B-1, (12) 


we see that the slopes of the optimal lines are —p/q and ~—gq/p. 
Thus, the equilibrium is stable if and only if p>gq. The case p = 
g = 1/2 was treated in I. There the optimal lines coincided, and 
the equilibrium was “‘neutral.”’ 

If each player is allowed to adjust his effort while the other 
keeps his constant, the behavior of the players will be as shown in 
Figure 1. Even if the adjustments are made simultaneously, the 
point (x, y) will move along the resultant of both adjustments, pro- 
vided no learning takes place. This means that no bias toward ad- 
justing the efforts upward is introduced by the fact that sometimes 
dz and dy are both positive and jointly bring @ and y to a position 
better for both (see below). The effects of learning would be in- 
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(a) (b) 

FIGURE 1. If the equilibrium is stable (a), the efforts of both players 
are so adjusted as to bring them toward the equilibrium. If the equilib- 
rium is unstable (b), the adjustment of efforts carries them away from the 
equilibrium until one of the players becomes a parasite. 


teresting to consider, but we will not pursue this matter here. We 
summarize our results in 

Theorem 2. /f the social matrix is diagonal and p > q, the ef- 
forts will be stabilized at x=y=p/B-1. If p <q, the equilid- 
rium is unstable. In this case, one player will become parasitic on 
the other; namely, tf the player to choose first takes a value less 
than p/B -1, he will become the parasite; otherwise the other 
player becomes parasitic. 

We note that neither the stable equilibrium satisfactions nor even 
always that of the parasite in the unstable case are the greatest 
that can be achieved. Indeed, if the pay-off function is known and 
the players enter a coalition, they can achieve the greatest total 
pay-off by maximizing S$, + Sg. In this case their individual S; will 
be greater than in the equilibrium. 


] P q 

Sas a8 ee ee eee = 5K 13 

je 2) Ll+pe+qy l1+qr+py B : (18) 

d q P 

so Spits, eames FE ci. mh st 14 

ay 2) L+pet+qy 1+qr+py Pp : ce 
which leads to the solution 


a* = y* = 1/8 -1. (15) 


Comparing now the resulting primary satisfactions with those pre- 
viously obtained, we have, substituting (15) into (5) and (6), 


Sf =S* =-logB-1+86, (16) 
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whereas (9) introduced into (5) and (6) gives 


S, =S, =logp-logB-p+B. (17) 
It is easily seen that for p< 1, S, < SF, since 
S, -S*¥ =logp-p+1<0. (18) 


In the parasitic case, if player I is the parasite, the ‘‘host” 
produces 


y=1/8 -1/p, (19) 


and the parasite’s satisfaction obtained by substituting (19) into 
(5) with 2 =0 is 


S** = log (1 + g/B - q/p) , (20) 


which may be greater or less than S*, depending on the values of p 
and 8B. Even if S** > S*, this gain occurs at a greater expense to 
the host so that the solution is not a “‘social optimum,”’ where S, + 
S, is maximized. 

The social optimum solution S* is given as a solution of our 
game in the classical theory, where the pay-off functions are as- 
sumed known to both players and where collusion is allowed. Itis 
interesting to note, however, that this solution can also be arrived 
at without knowledge and without collusion as an equilibrium, pro- 
vided the secondary satisfactions S$’ are introduced. (This was 
shown already by Rashevsky.) Indeed, the solution is clearly ar- 
rived at for all social matrices with A=B, C =D. It is also easy 
to show that those are the only matrices which lead to the social 
optimum solution. For setting the partial derivatives of (3) with 
respect to z and y respectively equal to zero, we obtain 


B 
ae ee ee ee (21) 
1 + pe + gy 1+ gx + py 
D 
A. EN clad Dre eT (22) 


Lliprt+gqy 1+qu+py 
If now the social optimum (1/8 — 1, 1/B — 1) is to satisfy (21) and 
(22), we obtain after substitution and re-arrangements 


B(A p- A+B g)=0=(B- A)BQ (23) 
B(C g-D+D p)=0=(C-D)Bq. (24) 


Unless g=0, we must therefore have A = B, C = D. The case g = 
0 also leads to the social optimum, but in that case the individuals 
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are totally independent of one another. As has been pointed out in 
I, some additional inducement is required for the establishment of 
mutual dependence. We take the existence of such inducements 
for granted so that the arrangement persists. 

In what follows, it will be useful to compute the equilibrium 
(z, 7) in terms of our parameters, A, B, C, D, p, g, and B. This 
can be done by solving the system (21) and (22) for ze and y. Itis 
simplest to solve first for the amounts received by each individual. 
We denote these by G/B and H/B, and so obtain 


G = ees = B(1 + pz + 99) (25) 
H = eect = B(1 + g@ + py), (26) 
from which we get 
E> gpg OG at -4 vat 
To genes (pH~4@)-1. (28) 


Equations (27) and (28) constitute the parametric equations of 
the path of the equilibrium as the social matrix changes. We note 
further that a necessary and sufficient condition that the path move 
along the line y = x is AC = BD, as can be verified by equating the 
right sides of (27) and (28) and substituting the values of G and H 
from (25) and (26). At the social optimum G = H = 1. 

Let now the two players begin with a diagonal matrix and let the 
process of “‘socialization’’ set in by having the elements B and C 
of the social matrix increase from zero to some positive value. We 
will confine ourselyes to the cooperative case, where all the ele- 
ments of the matrix are positive. Suppose also that p> q so that 
the initial position of (%, 7) is actually at the equilibrium (p/B — 1, 
p/B -—1). We are interested in the direction taken by the equilib- 
rium from its initial position as all the four elements of the matrix 
begin to vary. To do this we compute 


dy Sa Vee ea 29 
= ---- } ae GBS == — 
aA Bias edie k Oana uae ee = 
3 3 3 F , 
dy =. GA Ran GE nee og (30) 
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At the initial equilibrium, however, 
02/0A = dx/dD = dy/0A = dy/AD = 0. (31) 
This is a consequence of the invariance of the equilibrium with re- 
spect to the units of satisfaction when the satisfactions are inde- 
pendent. This can also be directly verified by differentiating (27) 


and (28) with respect to A and D partially, while setting B = C =0. 
Therefore equations (29) and (30) reduce to 


s Oz aB Ox dc 39 
fee dea te ee? o 
dy dy 
dy = --"- dB + ---- ; 
ee Olt 0 (33) 
We find these differentials to be at the initial point 
Gy aan 235% Yo) (34) 
BP - 9) 
1 
OY = ge (-qg?A dB +pqD" dC). (35) 
B(p — 9) 


From these equations we immediately derive 

Theorem 3. Let the initial equilibrium be stable. Then, as the 
cross terms of the social matrix acquire small positive increments, 
and regardless of the increments of the diagonal terms, player I 
will increase his effort <f and only if 


pegiA dB > g2D 7d, (36) 
and II will do likewise if and only if 
Pega YG As dB. (37) 


We now ask under what conditions will the respective primary 
satisfactions increase. Since at initial equilibrium 


AS,/de = OS8,/dy =0, (38) 
we have 

as | 
dS, = -=*. dy = Bq p™ dy , (39) 

oy 

aS, 
dS, = ae de=Bqp dz. (40) 

we 


We see that the primary satisfactions of each player will increase 
if the other raises his effort. The necessary and sufficient condi- 
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tion for both to enjoy a gain in primary satisfaction as the social 
matrix begins to depart from the diagonal is 


dx dy>0; dx+dy>O0. (41) 


From (41), (34), and (35) we obtain by pupae and yogis the 
right sides, and remembering that p + g = 1; hence p2+q?=1- 279, 


q° 


dx dy = eae [dB dC/AD - pq(dB/A+dC/DY] (42) 


dz + dy = g8'(dB/A + dC/D) , (43) 


and we have proved 
Theorem 4. /f q #40, a necessary and sufficient condition that 
socialization raises the primary satisfactions of both players is 


dBdC/AD > (dB/A + dC/DY pq (44) 
dB>0; dC>0. (45) 


If g=0, the equilibrium is already at the social optimum and so- 
cialization has no effect. The case p = q leads to indeterminacies 
and will be treated separately. 

Note: Although in this case an increase in the efforts of both in- 
dividuals was sufficient for the increase in primary satisfactions, 
this need not be the case when the equilibrium is unstable. In that 
case the initial situation is parasitic. It remains parasitic as long 
as the equilibrium remains unstable. Thus it is possible for the 
social matrix to change in such a way that the host increases his 
effort while the parasite remains inactive. But this forces the host 
off the line optimal with respect to his primary satisfaction. Thus 
his satisfaction can only decrease. Another example when an in- 
crease of both efforts may bring a decrease in both primary satis- 
factions is when the social matrix carries the (stable) equilibrium 
beyond the social optimum (as we shall see below). This decreases 
S, + S, as the efforts increase. If the equilibrium moves along the 
line y = x, S, = §,, each of the primary satisfactions decreases. 

Next we consider the effect of the social matrix on the stability 
of the equilibrium. The slopes of the optimal lines at their point 
of intersection are defined by the following expressions, obtained 
by differentiating (21) and (22) with respect to 2 and using equa- 
tions (25) and (26) 
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Now if the equilibrium is to change from a stable one to an un- 
stable one or vice versa, it must pass through a neutral state where 
¥; = ¥- Imposing this equality, we obtain after re-arrangements 
the following necessary condition for neutral equilibrium: 


ADp? = BCq?. (48) 
Referring, however, to the expressions for % and 7 given by (27), 
(28), (25), and (26), we note that if (48) holds, and 7 will be 
positive only if at least one of the denominators of G or H van- 


ishes. But (48) implies that if one vanishes, the other must too. 
Hence, in addition to (48), we must impose also 

Ap=Bq; Dp=Cq. (49) 
When (49) obtains, the optimal lines coincide and the equilibrium 
becomes indeterminate. It makes sense, however, to talk of the 
limits of % and ¥ as the optimal curves approach coincidence. This 
limit depends on the way Ap approaches Bg and Dp approaches Cq. 
They can, of course, do so independently. Suppose 


Bq=Ap+0,; Cq=Dp+6,, (50) 
where 0, and 6, approach zero. Then it is easy to show that 
Kp - K"q | 
Pig = h 5 a ee z | (51) 
6. + 0 B p-@q 
6, > 0 
K1p-K 
eae h ets Wes | (52) 
6,>0 B ae. 
6,0 
where K = Lim A@,/D0,. Since this limit is arbitrary if 6, and 6, 
yO 
0, > 0 


are independent, an arbitrary point can be reached when equilib- 
rium changes status. In the special case where the equilibrium al- 
ways stays on the line y =, we must have @ = 6,D/A, and the 
change of stability will occur at @ = y = 2p/B -1. 

Theorem 5. The equilibrium at the social optimum is always 
stable, except when p = q, in which case tt is neutral. 
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Proof. We recall that at social optimum, we must have A= 3, 
C =D. Imposing these conditions on (46) and (47), we obtain 


TA = ee (53) 


But unless p= q, we always have pg<1/4, therefore, 7 <-—1, 
and | ¥|>|%|-. Thus Theorem 1 implies Theorem 5. 

We can see the significance of this conclusion by noting that 
along the line y= 2, stability changes at r=Y=2p/B-1. If 
p> q, then p> 1/2, and the change takes place beyond the social 
optimum; hence, stability will obtain all the way. Contrariwise, 
the change will occur short of the social optimum as @ and 7¥ in- 
crease from their initial position of instability so that again the 
social optimum is stable. 

We now examine the case where p = gq. Certain aspects of this 
case were discussed in I. In particular it was argued that this 
case leads to parasitism, however, under the assumption that the 
first player to abandon the position p/B8 — 1 by lowering his effort 
‘‘learns’’ that as the other player adjusts his own effort, this act 
leads to a greater satisfaction for himself and, thus, he repeats it 
until his effort becomes zero. If learning is not allowed, the posi- 
tion of the equilibrium is indeterminate and remains so. Suppose, 
however, the social matrix becomes non-diagonal, i.e., socializa- 
tion occurs. The optimal lines now remain straight lines, namely, 


A+B=2 AB(1 + 2/2 + y/2) (54) 
C+D =2 DB(1 + 2/2 + y/2). (55) 


Unless AC = BD, the lines are distinct, and a glance at Figure 2 

shows that the player with the optimal line having the smaller 

y-intercept will become the parasite; in other words, player I will 
I 


0 


FIGURE 2. If p= q, the optimal lines are parallel. Thus in this case, 
socialization generally results in parasitism. 
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become the parasite if BD < AC. This is intuitively evident quali- 
tatively since BD represents the product of the ‘‘concerns’’ for 
player II and AC for player I. 

The formalism developed here has a certain resemblance to the 
abstractions of behavioral science, even as game theory does. The 
parameters p and g determine a miniature ‘‘economic system’? since 
they deal with quantities exchanged. The satisfaction functions 
S, and S, describe a miniature ‘‘psychology’’ since they determine 
the motivations of the individual. This psychology pertains, how- 
ever, strictly to the ‘“‘isolated’’ individual, the classical homo 
economicus, who, while participating in the social activity of com- 
modity exchange, is not influenced by the effects of his actions on 
the satisfactions of others. The social matrix imposes a social 
psychology on our system, where the individual is motivated not 
only by his own primary satisfactions but also by those of others. 
The primary satisfactions, however, remain distinguishable, and 
one can compute their variation as the secondary (social) satisfac- 
tions are modified by the social matrix. 

The model lends itself naturally to generalizations, for example, 
involving more than two individuals or more complicated systems 
of exchange, or else increased efficiency of output as a function of 
greater satisfactions, etc. Experimental procedures also suggest 
themselves. Human subjects could conceivably engage in a game, 
such as was described, where the primary and secondary pay-offs 
could be built into the rules without the rules themselves being 
made aware to the subjects. It would be interesting to see if, 
guided by trial and error, individuals would exhibit the effects en- 
visaged in the theory. One must keep in mind, however, that learn- 
ing is bound to occur with human subjects so that they may forgo 
immediate maximization to increase expected pay-offs in the future. 
A most interesting question is whether in the absence of verbal 
communication or even of contact of any kind, individuals can learn 
to cooperate so as to reach more advantageous points, even if these 
are not among the equilibria.* Or, one might make known to both 
individuals the primary satisfactions achieved by both, without, 
however, introducing secondary satisfactions into pay-offs. Will, 
under these conditions, cooperative behavior occur simply in view 


*Some pilot studies have shown this to be the case (personal communi- 
cation from H. Raiffa). 
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of the social habits which the individuals bring with them into the 
laboratory? 

Biological applications, although strongly suggested, would be 
quite difficult to reproduce experimentally because of the compli- 
cated nature of experimental materials. (We are not referring to be- 
havior experiments with animals, in which there has developed an 
extensive technique, but experiments involving metabolically 
founded symbiosis and parasitism.) Nor is it expected that a reali- 
zation of so simple a model could be found in the biological world. 
The value of the model is confined to the conceptualizations sug- 
gested by it. 


I am indebted to Professor Howard Raiffa for his critical and 
helpful discussion of this paper. 
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The twentieth century has witnessed a geometrization of physics, that 
is, a reduction of the basic concepts of physics to geometric concepts. 
The topological approach to biology, recently proposed and to some ex- 
tent developed by the author, is a small step in the direction of geo- 
metrization of biology, but is unable to achieve the main purpose of such 
a geometrization of biology, namely, the reduction to geometric concepts 
of such purely biological concepts as ingestion, digestion, assimilation, 
etc. To achieve this purpose we must find geometric structures or spaces, 
in which different geometric properties stand to each other in the same 
formal logical relation, as the different concepts of biology stand to each 
other. If this were possible, then a set of geometric theorems could be 
‘‘translated’’? by an appropriate ‘‘glossary’’ into a set of biological laws. 

While not offering a solution to this problem, the present paper illus- 
trates the possibility of such an approach on several examples. Certain 
new types of topological spaces are introduced, which are used for il- 
lustration purposes only. It is shown, however, how from a theorem about 
such spaces a verifiable biological prediction could be made, if these 
spaces were to be taken seriously. 

A possible application to biology of E. Artin’s theory of braids is 
outlined. 


The study of mechanics began much earlier than that of any other 
branch of physics. Although Newton was the father of contemporary 
analytical mechanics, yet different mechanical laws, like those of 
the levers or those of flotation of bodies, have been known and 
were studied since the days of Archimedes. This historical situa- 
tion developed among physicists the attitude to consider mechanics 
as the basic branch of physics. The fundamental concepts, and 
even the fundamental laws of mechanics, became to be regarded as 
something self-evident. Newton referred to his laws as the *‘axioms 
or laws of motion,’’ and it was only at the turn of this century that 
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Henri Poincare pointed out that these were neither self-evident 
axioms nor experimental laws, but what he called useful conventions. 

Of the other branches of physics, the phenomena of sound were 
relatively early and easily shown to be of a mechanical nature. 
Thus acoustics became a branch of mechanics. The second half 
of the nineteenth century witnessed a similar ‘‘mechanization”’ of 
thermodynamics. Attempts to ‘‘explain’’ in mechanical terms the 
phenomena of light and of electromagnetism, especially the latter, 
proved, however, to be unsuccessful. When it was demonstrated 
that propagation of light is a periodic phenomenon, it was natural 
to interpret it as a wave motion in some invisible and even “‘im- 
ponderable’’? medium, and the concept of universal ether was cre- 
ated. As is well known, this concept became more and more fraught 
with difficulties and self-contradictions. 

Yet all the great physicists of the 19th century believed—nay, 
were almost convinced—that all phenomena of physics are me- 
chanical in their nature. Faraday thought of the electromagnetic 
field in mechanical terms; Maxwell found his equations by a series 
of unsuccessful attempts at interpreting electromagnetic phenomena 
as manifestations of mechanical ones. The equations were estab- 
lished and proved their worth by their enormous predictive value 
only, however, to show that they were the inspirational guess of a 
genius rather than a logical consequence of a mechanistic picture, 
Yet, outstanding physicists like L. Boltzmann and others continued 
the unsuccessful efforts to reduce Maxwell’s equations to Newton’s 
laws of motion. 

Inasmuch as Maxwell’s work established the identity of optical 
and electromagnetic phenomena, showing the former to be only a 
branch of the latter, the difficulty of interpreting mechanically 
Maxwell’s equations only further increased the difficulties of a me- 
chanical interpretation of optical phenomena. 

Towards the turn of the century the notion began to crystallize 
that perhaps electromagnetic phenomena not only cannot be reduced 
to mechanical ones, but that actually they are the basic phenomena 
of physics to which all other phenomena, including mechanics, may 
be reduced. This idea was especially suggested by the discovery 
of the electromagnetic mass of the electron, predicted theoretically 
by H. A. Lorentz and subsequently verified by several experimenters. 

The notion, however, did not look to be a happy one to the 19th 
century classical physicists. The notions of electromagnetic field, 
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electric and magnetic vectors did not have the property of immediate 
intuitive self-evidence, which the mechanical notions of motion, 
mass, acceleration, and force seem to possess. It was not under- 
stood that this self-evidence of the mechanical concepts was only 
an apparent one, due to the centuries of familiarity of the physicists 
with those concepts. This self-evidence was nothing more than a 
habit of thinking; but even the great scientists have frequently 
found it difficult to free themselves from long established habits of 
thinking. The laws of human psychology are very much the same 
for the scientist and for the layman! 

It was in 1909 that Hermann Minkowski, inspired by the profound 
discovery by Hinstein in 1905 of the relativity of space and time, 
indicated a way of reducing physical phenomena to something 
which is outside of physics, even outside of any other natural 
science, namely, pure geometry. Minkowski showed that the basic 
concept of mechanics, and possibly of all physics, the concept of 
motion, may be interpreted as an orthogonal transformation of a 
system of coordinates in a four-dimensional hyperspace with one 
imaginary coordinate. Inspired by Einstein, this discovery in its 
turn eventually led Einstein to the creation of the General Theory 
of Relativity, in which not only the concept of motion, but also’ 
concepts of acceleration, mass, and force, and especially the no- 
tion of gravitation received a geometrical interpretation in a four- 
dimensional non-Kuclidean hyperspace. 

True enough, a complete geometrization of physics is still a 
relatively distant goal. One of the stumbling stones is again Max- 
well’s equations. However, the difficulties now appear to be of a 
different kind; they are more of a mathematical than of a conceptual 
nature. Beginning with the ingenious—though unsuccessful at- 
tempt by Herman Weyl (1920)—attempts at a reduction of electro- 
magnetic phenomena to geometric concepts have continued, and 
Einstein’s own later contribution (Einstein, 1953) seems to carry 
a great deal of promise. 

But even to the limited extent that geometrization of physics has 
been achieved, it proved to be important. Not only has it shown an 
unusual predictive value both in physics and astronomy, but it is 
different from the old attempts to reduce physics to one of its 
branches in that the reduction of actual natural phenomena is made 
to purely mathematical concepts. Those concepts are creations of 
the human mind, and basically are therefore much more intuitive 
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than the physical ones, even though a visualization of multidi- 
mensional spaces may be barred to us. Basically, what modern 
physics does, is to map the observable physical phenomena iso- 
morphically onto an abstract geometric structure. Certain concepts 
of physics, like motion, acceleration, etc., map on purely geo- 
metrical concepts of coordinates, curvature, etc. The theorems of 
geometry, which establish metric relations between different geo- 
metrical concepts, then lead directly to laws of physics by the use 
of .a sort of glossary or dictionary which shows the name of the 
physical concept that is to be substituted for the name of a geo- 
metric concept. A ‘‘translation’’ of physical laws into geometric 
theorems and vice versa is thus possible. 

Compared with physics in age, biology is almost a newborn baby. 
Whenever we experience outside of ourselves any manifestations of 
life we do it only through its physical manifestations. Certain 
introspective psychological experiences, as well as the more im- 
mediate contact with life than with the nonliving, may have led 
some biologists to assume that phenomena of biology are basically 
different from those of the nonliving physical world. An utterly 
useless and time-wasting argument ensued between the vitalists 
and the mechanists, which are now better called physicalists, as 
to whether phenomena of life are something sui generis, nonphysi- 
cal in nature, or are basically reducible to the laws of physics. As 
remarked above, the only objectively scientific study of life can 
be made through the study of its physical manifestations. This 
holds also for the so-called mental phenomena. We can know what 
another person thinks only by his telling us, or writing it, which 
are both physical acts, or we may infer about his thoughts from 
some other of his overt behavior, which again must be physically 
manifested if we wish to observe it. Therefore the scientific study 
of biology has’ become a study of different very special physico- 
chemical situations, with the applications of the methods of both 
experimental and theoretical physics. This approach has proved 
to be of tremendous success. As to whether all phenomena of 
biology can be explained in terms of contemporary (1955) physics 
is a question which cannot be answered by any general specula- 
tion. All we can say that attempts at such explanation, when made 
competently, have been hitherto crowned with success. It is, how- 
ever, impossible to deny that a biological phenomenon may be dis- 
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covered which cannot be explained in terms of physical laws known 
at present. But far from proving anything ‘‘unphysical’’ about life, 
such a discovery will merely mean the need of an extension of 
physics (Rashevsky, 1934, 1955b), just as the impossibility of ex- 
plaining some spectroscopic phenomena in terms of classical phys- 
ics has led to its extension by introducing quantum mechanics. No 
one will ever call quantum mechanics ‘‘unphysical.”’ 

If physics, in its present or in its extended form, is to be re- 
duced to geometry, then it follows from all the above that eventually 
a geometrization will be also the fate of biology. When we come 
to consider this possibility more closely, we notice one very es- 
sential difference between some laws of biology and the laws of 
physics. All phenomena of physics are quantitative in nature. 
They not only can be measured, but even more, unless they are 
measured, very little if anything significant can be said about them. 
Many biological phenomena are also quantitative in nature, and the 
biologist is now becoming used to both measurement and mathe- 
matical analysis. But very many biological phenomena, and per- 
haps the most basic of them, are not quantitative but relational. 
Yet very definite statements can be made about them, and their 
importance is unquestionable. 

As has been pointed out elsewhere (Rashevsky, 1954), the re- 
sponse to a stimulus by a paramecium, followed by a locomotion 
toward food; the ingestion of the food, followed by digestion; the 
resulting absorption and assimilation of the digested products on 
one hand, and the defecation of undigested residue on the other— 
all those individual biological processes stand in basically the 
same relation in the paramecium as they do in a human being. The 
processes of stimulation, locomotion, digestion, etc., in a higher 
organism are much more complicated than in a protozoan and con- 
sist of many more partial ‘‘subprocesses.’’ But the corresponding 
complicated processes in a higher organism can be mapped in a 
many-to-one way onto the processes of a one-celled organism in 
such a manner that the basic relations remain invariant. 

If we wish to describe mathematically the above situation, we 
need a different mathematical apparatus from the one used hitherto 
in physics or in mathematical biology. Such an apparatus is pro- 
vided by topology, which is a geometry of relations rather than of 
quantities. To use this apparatus we must find some method of 
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describing an organism in terms of some topological spaces or 
complexes. The word ‘‘organism’’ itself suggests that a graph 
used as an organization chart of the organism may be employed to 
advantage. Such considerations led us to the formulation of a 
general principle, the principle of bio-topological mapping (Rashev-- 
sky 1954): 

The graphs which represent the organization charts of dif- 
ferent organisms are such that they all can be mapped on the 
graph or graphs of one, or at most very few, primordial or- 
ganisms. The graphs of all organisms are obtained from the 
graph of the primordial organism by the same geometrical 
multiparametric transformation, the graphs of different organ- 
isms corresponding to the different choices of the parameters. 


It has been shown that a topological biology developed on that 
basis with the addition of some special hypotheses has a definite 
predictive value. Not only is it possible to derive the known facts 
—that the more complex organisms possess a lesser regenerating 
ability for lost organs (Rashevsky, 1955c) and possess a greater 
adaptability (Rashevsky, 1955d), but it is also possible to derive 
the total number of possible different organisms (Rashevsky, 1954; 
1955a). The number thus found is about 100 times the actual num- 
ber of known species. Whether this means that nature has pro- 
duced hitherto only about 1% of all possible organisms or that 
some of the special hypotheses must be modified within the gen- 
eral principle cannot be said at present. The important thing, 
however, is that topological biology predicts something that can 
be verified by observation and thus offers a stimulus to new ex- 
perimental work. This is even better illustrated by the prediction 
of a relation between the number of tissues, the number of organs 
and the number of distinct biological functions (Rashevsky, 1955d), 
as well as by an estimation of the total number of hormones to be 
discovered in organisms of different degree of differentiation (Ra- 
shevsky, unpublished). An extension of combinatorial topological 
considerations to organic molecules (Rashevsky, 1955e; Karreman, 
1955) has led to biochemical conclusions which are in principle 
verifiable. 

The above approach is to some extent a move in the direction 
of a geometrization of biology. How does it relate to the geometri- 
zation of physics? Topological considerations have hitherto 
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played no part in the discussion of the four-dimensional hyper- 
Space onto which the physical phenomena are isomorphically 
mapped. The metric properties of the four-dimensional hyperspace 
is all that mattered. In many instances the metric of a space de- 
termines its topology. If this were true generally, we would expect 
that topological relations of the four-dimensional universe would 
be derived from its metric, such that they may lead to the descrip- 
tion modo geometrico of the basic biological relations. We then 
would have the following possibility of reduction: biology —> phys- 
icS —> geometry. However, as we have seen on a previous example 
(Rashevsky, 1954, p. 345), a situation may be conceivable in which 
life is a manifestation of ‘‘local’’ topology of the universe, in such 
a way, that while biological phenomena will always follow the 
laws of physics, they may not be reducible to those laws, while 
both physics and biology will be reducible to geometry. 

However, the geometrization of biology, as attempted by topo- 
logical biology, is still something very different from the geometri- 
zation of physics. In the former we do not interpret in geometrical 
terms any of the biological functions, such as digestion, assimila- 
tion, etc. We may perhaps say that we eliminate from biology the 
elusive concept of ‘‘organization,’’ and substitute for it the geo- 
metrical notion of the topological properties of the graph of an 
organism. But except for that, we still manipulate with such con- 
cepts as digestion, assimilation, locomotion, stimulation, etc. We 
consider those concepts as given, and all we do is to study topo- 
logical relations between those concepts, in an abstract space. 

The geometrization of physics goes much further. It interprets 

‘the basic concepts of physics in geometrical terms. It eliminates, 

in a sense, those concepts from physics. It does not merely seek 
formal geometric relations between different masses, velocities 
and accelerations. This has been done long ago by classical 
kinematics and does not constitute a mapping of physics onto 
geometry. 

If we wish to achieve the same thing for biology, we must go 
much further than topological biology has done so far. It seems 
that combinatorial topology is hardly the adequate apparatus for 
our purposes and that recourse to set topology will have to be 
made (Rashevsky, 1955d). 

The whole universe, physical and biological, is a set of “‘ele- 
ments,’’ the word element not being used in the chemical sense. 
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In order to avoid confusion with the chemical connotation of the 
word element, we shall use it in quotation marks when set theoreti- 
cal connotation is implied, and without quotation marks when the 
chemical connotation is used. Perhaps the most basic way of 
looking at it is to consider the universe as a set of world lines. 
We may, however, consider as ‘‘elements’’ of the sets, groups of 
the world lines, or their points of intersections, or other configura- 
tions. Different chemical reactions may be considered as the 
‘“elements’’ of the set. In any case, however we define the ‘‘ele- 
ments’? of this set, the biological phenomena are a subset of it. 

What are the characteristic properties which make us recognize 
an organism as such? Perhaps the most basic thing about an or- 
ganism is that it selects certain ‘‘elements’’ from the inorganic 
environment even if those elements form parts of other subsets. 
Once a selection is made, those selected ‘‘elements’’ are organ- 
ized in a definite pattern eventually becoming themselves ‘‘ele- 
ments’’ of the organism, and leading to a duplication of the latter. 
In this case the ‘‘elements”’ are the different chemical molecules, 
radicals or complexes. 

This process of selection is accompanied by a loss of some 
‘‘elements’’ of the organism. This loss may be a complete break- 
down into waste products, or it may be constituted by the secretion 
of a molecule of a digestive enzyme which ‘‘attacks’’ a food par- 
ticle and selectively breaks off it the desired molecular configura- 
tion. Usually the secretion of a digestive enzyme is not con- 
sidered on the same level as a catabolic breakdown. But logically 
they both resemble each other, and they both may be necessary for 
the process of selection. For while in some microorganisms, like 
nonmotile bacteria, the selection goes on at a purely chemical 
level; in higher organisms it is more complicated. The sight of 
possible food; the locomotion towards it; the breaking up of food 
either by teeth, or by hands, or by hand-made machinery; the pick- 
ing up and consumption of the proper parts of the food—all this 
constitutes a process of selection, and for many such processes 
in higher organisms the energy released by catabolic processes is 
needed. Thus everywhere in the process of selection some loss 
of its constituents by the organisms seems essential. 

Any selection process implies rejection of the unwanted material. 
Digestion, followed on one hand by absorption, on the other by 
defecation are only the logical aspects of the selection process. 
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Thus these three very basic biological functions may be described 
in terms of selection. And as we have just said, a number of other 
biological functions, which acquired typically biological names, 
are also basically logical aspects of the process of selection. 

If we could interpret selection geometrically—or more specifi- 
cally topologically—we would thus reduce some biological con- 
cepts to geometric ones. 

Let us make the next logical step. Selection implies the divi- 
sion of the ‘‘elements’’ of some subset in two classes: those se- 
lected and those rejected. Division in two classes is not at all 
uncommon in topology. Thus a point of the one dimensional space 
of real numbers induces a Dedekind cut, and, therefore, a division 
of all other points in two classes. A closed Jordan curve in F ? 
divides all points of £? in two classes: the inner and the outer 
ones. Both cases are examples of a subset ¥ of a space S dividing 
S~—M in two classes. But those simple cases do not help us any. 
As we have seen, the division by the organisms of the set of ‘fele- 
ments’’ which constitute the environment of an organism in two 
classes is contingent on the organism, considered as a set, losing 
some of its own ‘‘elements.’’ Moreover, if one or more of the se- 
lected ‘‘elements’’ are added to the subset M which is the organ- 
ism, the relation between that subset and the complement S — M re- 
mains the same: the organism continues to select. That second 
property is exhibited in the following example. Let M be a subset 
of S not closed in S. Then M divides all the points of S — M in two 
classes: those that are the limit points of M, and those that are 
not. Except for the very special case in which the first class is 
a degenerate subset, when some of the limit points of / are added 
to M the resulting subset M’ still divides the points of S— M’ in 
two classes: those that are limit points of M’, and those that are 
not. 

A particular case of the above offers some remote analogy to the 
first necessary property, that of a loss by breakdown. Let S be 
the space of real numbers and let M be the subspace obtained from 
S by omitting all rational points in the closed interval (0,1). All 
rational points of S within that interval are limit points of M; those 
outside are not. The subspace M is not connected, however, and 
this loss of connectivity may be regarded as a remote analogy of 
a breakdown process. In this example there is no division by MM of 
S —M in two classes, that is, no selection, without loss of con- 


nectivity of M. 
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The above examples have, of course, no scientific value what- 
soever. They are used only to illustrate a point which, if given 
an abstract formulation, may be difficult to make clear. Namely, 
if we can find a topological space, such that its different topo- 
logical properties stand in the same logical relation to each other 
as do different basic biological concepts, then by substituting the 
names of those biological concepts for the names of the corre- 
sponding topological properties, we could, using again a glossary, 
translate every theorem about such a space into a corresponding 
biological law, just as it is done in the geometrization of physics. 
The purely mathematical study of the properties of such a properly 
chosen space would yield us all the laws of biology, both already 
known and those not yet discovered. If this ever becomes pos- 
sible, it would indeed be a triumph of geometry, more generally, of 
pure mathematics. 

This idea, suggested by the twentieth century physics, may ap- 
peal to some and not appeal to others. But even to those to whom 
it may appeal, it still is of no use. You cannot do much with a 
bare idea unless you at least show how, in principle, the idea 
might be realized. Here we must consider two logical possibilities. 

Topology is still a young science; yet, it has already ramified 
in many directions, and a very great variety of its aspects have 
been studied. It is not precluded that the necessary types of 
spaces or other topological structures have already been studied 
and only need to be taken over, but that thus far we have failed to 
perceive the proper possibility. After all, the basic apparatus of 
non-Kuclidean geometry has existed for a fairly long time; yet, it 
took Einstein to see that phenomena of physics are isomorphic to 
some of those already well known geometric properties. The ra- 
tional problem-solving and the reduction of new situations to al- 
ready known ones has psychologically and biophysically some 
similarity with the discovery of a ‘‘hidden picture’’ as has been 
shown elsewhere (Rashevsky, 1948; Chapt. XLII). Thus, we may 
hope that an inspiration by a mathematician or a mathematical 
biologist may indicate the solution, which will appear rather ob- 
vious after it is discovered. 

But we cannot discount the other possibility, namely, that ap- 
propriate topological spaces have not yet been studied and that 
the needs of biology may provide a stimulus to the topologist to 
make new discoveries in his field, just as the needs of physics 
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have provided in the past the stimulus for new mathematical dis- 
coveries. Gauss has actually credited physics with many of his 
purely mathematical inspirations. 

But here again we must be somewhat more definite and therefore, 
even without hoping to formulate at present the necessary purely 
mathematical problem and to describe the necessary type of topo- 
logical space, give at least a few examples of the kind of thinking 
that may be required. Those examples we shall give again as an 
illustration only. Their actual biological and mathematical values 
are nil. But they may serve, as the French say, ‘‘pour fixer les 
idées.”’ 

Let us return again to the basic logical property of the organism. 
Certain ‘‘elements’’ of the environment are selected, that is, begin 
to belong to one of two classes if some ‘‘elements’’ of the organ- 
ism itself are detached from the latter. Considering again the or- 
ganism as a subspace M in the space S and considering the ‘‘ele- 
ments’’ as points in these spaces, we may translate the above 
logical property into topological language. In order that M should 
induce a division of the points of S — M in two classes, it is nec- 
essary that some points of § ~M would possess some property P 
with respect to M, which others do not possess. In order that the 
division of the points of S —M in two classes becomes possible 
only after M loses some ‘‘elements,’’ some points of S-— M must 
acquire the property P with respect to M only after certain points 
of M are removed. The simplest example of such a situation is 
offered by a closed Jordan curve, from which one point, « (Figure 
1), is removed. Let such a Jordan curve with a point removed be 
the subspace M in the space S=E?. It is actually a Jordan arc, 
the ends of which are separated by the point « only. Thus M isa 
single connected subspace of FE”. The point ~ does not belong to 
M, but belongs to S—M. Let property P consist of making M con- 
nected. Since M is already connected, therefore « does not pos- 
sess property P. But if we remove from M any other point, say B 
(Figure 1), then M becomes non-connected, and « acquires the 
property P. 

Instead of a closed Jordan curve with one point removed, con- 
sider the structure shown in Figure 2, in which the points o,, 4, 
and a, are removed. As a subspace M of S = E”, it is connected. 
If, however, we remove the point 8 from M, it breaks up into four 
components, and the set of points (4, Oy, X,) acquires the prop- 
erty P. Thus the removal of one point, B, makes M divide the 
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different subsets of S—/M in two classes: one class, consisting 
only of one element—the subset (a,, &,, 4,)—has the property Pe 
the other does not. The subsets of the two different ‘‘classes”’ 
are not all disjoined. Thus the subset (a,, d,, 4,) does not have 
property P, but (o,, O,, 4) > (ys Ags %,) # 0. 

In a similar space like the one shown in Figure 2, but having in- 
stead of 4 lines n + 1 lines connected by the point 6, the removal 
of the point B will make a set of n points of S—M acquire the prop- 
erty P, where n is an arbitrary positive integer. 

We shall call such subspaces of E”, as those shown in Figure 2, 
A,- spaces. Now connect an infinite number of A,|— spaces in a 
manner shown in Figure 3, in which all A, =A,. Considered as 
subspace M of FE? = S the space thus obtained is connected, non- 
compact, dense in itself and nowhere dense in S, and not closed in 
E?, since all the « points are its limit points. It satisfies the four 
Hausdorf axioms. It has the property that the removal of any B, 
makes a set of points (con oe ... a") acquire the property P. Any 
subspace of S that has this property we shall call a » - space. 
The number of points 8, is countably infinite, and so is the number 
of sets of S — M that can acquire the property P. 

We can make the A — space shown in Figure 3 bounded in S = EF? 
and derive a more complicated }— space. To achieve this, take 
the spiral with the equation in polar coordinates p, 0: 


Oe Pa 
Dg (1) 


At equal, or unequal intervals, along the line (1), which begins 
at p=0 and approaches asymptotically the circle p = 1, attach to 
it A — spaces, as they are attached in Figure 8 to the line AB, 
(Figure 4), making the size of each A - space smaller, as the 
windings of the spiral become closer, so that no A — space would 
be intersected by the line (1). Then take spirals given by: 


o-(t +m] 

el iA sia 3 0 > 1+ < (2) 
where m is an arbitrary integer different from zero. The spirals 
(2), obtained for different values of m, lie all between the windings 
of spiral (1) (Figure 5). Attach to aes spiral (2) a A, — space, 
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at finite intervals, decreasing the size of A, — spaces as the m 
increases and as the spirals come closer to each other. If m in- 
creases to infinity, in the limit, we obtain a A — space, which lies 
all within the circle of radius 1. It has a countably infinite number 
of 6 points and there is a countably infinite number of sets of 
points of S — M (M =); S = E?) which can acquire property P. The 
\ — space in question is not compact because a sequence of equi- 
distant points along any of the spiral lines does not have a limit 
point. As subspace M of E?, the d - space is connected, dense in 
itself and nowhere dense in FE? = S. It satisfies the four Hausdorf 


axioms. 
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In the \— space considered above, the cardinal number of the 
set of B — points, as well as of the set of sets that can acquire the 
property P, is %,, whereas the cardinal number of the points of the 
whole A ~ space is ¥, that of a continuum. 

We now shall construct a different space M, in which the re- 
moval of one point makes an infinite set of points of S ~ acquire 
property P. Consider the family of circles in E?: 


(2~o)? +y? =a, (3) 


where « = % or & is an irrational number between 0 and 4. 

The cardinal number of the set of such circles is 8. Each circle 
intersects the line y =0 at 2 =0 and at the point x= 2a. Remove 
in each circle the point x = 24. The space M obtained by this re- 
moval from the family of circles (3), considered as a subspace of 
S =E7, is connected. Removal of the point 0, which we shall 
designate by B, from M makes it not connected, and in order to 
make it connected with the point 0 removed, we must add to it all 
the points 2 =2a. Hence the set of all irrational points in the 
interval (0,1), whose cardinal number is 8, acquires now property 
P. We shall call this space M a X’ — space. 

Let us now consider the family of circles (3) not in the plane 
2 =0, but in a plane gz = irrational number, and 0 <z2<1. Remove 
in each plane all the points of intersection x = 2a and connect the 
outermost circles, those corresponding to “=%, by any line 1 
parallel to the z ~ axis, except the 2 — axis itself, passing through 
that outermost circle. 

Such a set of circles with points removed as indicated, con- 
sidered as a subspace M of E°, is connected and bounded in E?. 
The removal of any point 8 with irrational z ~ coordinate on the 
line 7 destroys the connectedness of M, and the set of points 2a 
which correspond to that zg-— coordinate acquires the property P. 
The space so constructed we shall call the A’’ ~ space. 

The following theorems can be readily proved about the X, \’ and 
d’’ — spaces. 

The intersection of two i, (A’, \’’) spaces is not a A, or 
\”” space, and is always not connected. 

As we have emphasized already, the above examples of ‘‘fancy”’ 
topological spaces are of no value either biologically or topologi- 
cally and are used only to illustrate a point. We shall now use 
them for such an illustration. 
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An organism, as we have seen, selects from the environment 
some “‘elements’’ which have some special properties with respect 
to that organism. The selection of a set of ‘“‘elements’’ becomes 
possible only after the organism loses some other set of its own 
‘“‘elements,’’ and, in general, a specific set must be lost in order 
to make the selection of a given set possible. The selection, as 
we also have seen, is logically equivalent to inducing the division 
of all the ‘‘elements’’ of the environment in two classes: those 
that possess certain properties, and those that do not. 

The A, A’, and A’’ — spaces are, from a logical point of view, se- 
lecting the specific sets from the environment § —M upon losing 
other specific sets. We have here a description of a biological 
phenomenon in topological terms. Or we may put it this’ way. If 
we make the ‘‘elements’’ of an organism, which must be lost in 
order to make a selection from the environment, correspond to the 
B — points of the X, »’, or A’’ + space and the selected ‘‘elements’’ 
to the «— points, then the logical relation between the loss of 
*‘elements’’ by an organism and its selection activities are mapped 
isomorphically on the topological relations between the aand B 
points. 

In an organism the loss of an ‘‘element’’ is compensated by the 
addition of the properly selected ‘‘elements’’ of the environment. 
If, after we remove a $8 — point, we add the corresponding %- 
points, we compensate by this addition for the loss of the B — 
point and restore the original connectedness of the A, 4’, or A’’ — 
space. 

Now, to continue the illustration, let us just for a moment imag- 
ine that the A, \’ and A’’ — spaces would be used not only for il- 
lustration, but so to say, taken seriously. 

The selection of ‘‘elements’’ of the environment through loss 
of some of its own ‘‘elements’’ is a basic property of the organism. 
We have seen that an intersection of two X, 4’, or \”” — spaces is 
not a connected space, and it does not possess the property of se- 
lection which the A, \’, and A’’ — spaces possess. Let us trans- 
_ late the theorem into biological language. It reads thus: the com- 
mon part of two organisms is not an organism and does not itself 
select the proper ‘‘elements’’ from the surroundings. Since a cell 
is an organism, we may restate the above as follows: The com- 
mon part of two cells does not possess the ability of selecting 
proper elements from the environment by losing some of its own 


elements. 


66é 
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But the only possible common part of two cells could be the 
cell membrane which separates two adjoining cells. Hence from 
topological considerations, we reach the biological conclusion that 
the cell membrane in such cases does not make a selection of 
proper material from the environment in a manner which the cell as 
a whole, and possibly some of its parts, do. Whether this is true 
or not, it is certainly an experimentally verifiable conclusion. 

The above example illustrates how the approach suggested in 
this paper, if properly made, can lead to verifiable predictions. 

Very many objections can be raised against the A, A’, and A’’ — 
spaces. Inasmuch as they are used only for illustration, we do not 
need to worry about the numerous possible objections. One point, 
however, may be mentioned perhaps as possible suggestion for the 
directions future research may take. 
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FIGURE 6. 


The number of ‘‘elements’’ of which an organism is composed is 
very large but finite, whereas the number of points in the A, )’, and 
\’’ — spaces is infinite. Moreover in the \ — space the cardinal 
number of all points is 8, while the cardinal number of 8 — and 
4—- points is ee We may say that only an infinitesimal fraction 
of the total number of points of the \ ~ space take part in the ‘‘se- 
lection’’ process. This is certainly not so ina biological system. 

It seems to be desirable to study spaces with a finite number of 
points. As an example used again for illustration purposes only, 
let S be the space formed by the points in a plane arranged in an 
infinite square lattice (Figure 6) and in which the neighborhoods 
are defined as follows: the neighborhood of a point consists of 
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that point and of the four adjacent points which lie on the lines 
of the lattice. Thus the neighborhood of point 1 in Figure 6 is 
constituted by the points 1, 2, 3, 4, 5. With neighborhoods so de- 
fined, S is connected because there exists no partition § = 4|B 
such that An B=AnNB=0. The subspace M, which consists of 
the points marked by heavy dots, is also connected and has a 
finite number of points. The removal of point £, from M makes the 
points a, and a, of S — M acquire the property P; while the removal 
of the point 8, makes the points a} and a3 of S~M acquire the 
property P. More complicated spaces of that kind can be readily 
constructed. 
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FIGURE 7. 


The following last example shows the possibility of having 
spaces which, as a result of a selection followed by addition of 
the selected ‘‘elements’’ (assimilation in biology), may duplicate 
themselves. The biological interest of such spaces is quite ob- 
vious. 

Consider in S = E? the space M shown in Figure 7a. From the 
point of view of combinatorial topology, M is a tree with the bi- 
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centers 6 and e (Konig 1936, p. 64). We shall consider it, however, 
from a set topological point of view as a subspace of S = E?, Mis 
irreducibly connected about its subspace 


C=AUB, (4) 
where A and B denote respectively the segments ac and df: 
A = ac; 'Brs dj. (5) 


Hence the segment be, which we shall denote by D, may be defined 
as 


De ee (6) 


The end points of D are cut points of A and B respectively. They 
divide A and B each in two components: A’, A’’, and B’, B’’. Ifwe 
remove from M either a point or a nondegenerate subset of D which 
does not contain the cut points of A and B, then M loses its con- 
nectedness, and any Jordan arc in § = E® which joins the sets A 
and B has the property P. In Figure 7b, line 6’e’ has, for example, 
the property P. 

Let us, however, remove from M two points 6, and 6,, which 
are not end points of M and which are contained in A’, A’’, B’, or 
B’’, and such that 6, and f, are not contained both in the same 
one of the above four subspaces. Thus in Figure 7c, B, is con- 
tained in A’’, 8, in B’; in Figure 7d, 8, is contained in A’, B, in 
B’. The removal of 6, and 8, from M separates M into three com- 
ponents, one of which contains D [expression (6)] and the other 
two which do not contain it. Denote the first component by M’, the 
other two by M’’ and M’’’. The component M’ is homeomorph with 
M; the components M’’ and M’’’ are not. 

Now, after the removal of 8, and 8, from M, any Jordan arc in 
S =E? which connects M’’ and M’’’ and whose end points are cut 
points of M’’ and M’’’ has the property of reconstituting a second 
space M,, which is homeomorph to M. After performing thus a 
specified operation, which involves the loss by M of two properly 
specified points 8, and 6, and the addition of a Jordan arc from a 
set of arcs selected by the loss of 8, and B,, we now have two 
spaces M’ and M,, each homeomorph to the original M. In Figures 
Tc and 7d such arcs are, for example, b’e’. 

Thus, removal of two properly located points from M imposes on 
certain Jordan arcs of § ~ M the property of duplicating the space 
M. The process can be repeated again with M’ and M, indefinitely. 
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It is true that the sizes of the subspaces which correspond to A’, 
A’’, B’, and B’’ decrease with each ‘‘reproduction.’? However, we 
are interested here not in the metric, but only in the topological 
properties. And topologically all the ‘‘daughter spaces’’ are ho- 
meomorph to the original space. Moreover, in spite of the reduc- 
tion in size, the cardinal number of the points contained in the 
‘daughter spaces’’ is always X. 

By study of appropriate topological spaces, it thus may be pos- 
sible to map the properties of the simplest conceivable organism 
onto a proper topological space. Using then the principle of bio- 
topological mapping, it may be possible to construct, by proper 
transformations, more complicated spaces, which map continuously 
in a many-to-one manner on the simpler space, and the study of the 
topological properties of those spaces may lead us to the discovery 
of new biological properties of multicellular organisms. 

Continuing our fantastic excursion into the possible future of 
topological biology, we shall now take a look at another branch 
of topology which, perhaps, may also be destined to contribute its 
share to mathematical biology. 

With the present-day reduction of physics to geometry, the phys- 
ical events are represented geometrically by the intersections of 
the world lines of different particles. The metric characteristics 
of those intersections, or as the physicists call them, space-time 
coincidences, describe also the metric characteristics of the phys- 
ical events. Inasmuch as an organism is composed of physical 
particles and obeys the laws of physics, the living phenomena are 
also representable by the intersections of world lines of the par- 
ticles of which the organism is built. However, the characteristic 
and basic properties of life being of a relational character rather 
than of a metric one, it is the topological relations between the 
intersections of world lines that are important now. The branch of 
topology which comes close to the study of such types of relations 
is the theory of knots (Reidemeister, 1932) and especially its sub- 
division, E. Artin’s theory of braids (Artin, 1925). 

Artin studies braids in a three-dimensional space, but an ex- 
tension to four dimensions seems to be natural and worth trying. 
Since again we shall use here some elementary notions of the 
theory of braids only for purposes of illustration, we shall confine 
ourselves to the three-dimensional case, studied by Artin. In fact, 
we shall consider even for Sep city only a degenerate two- 
dimensional case. 
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A braid is basically a set of lines interwoven in a regular man- 
ner. Consider the three lines of Figure 8. Let them first run 
parallel and then, at a certain point, begin an orderly interweaving. 
Such an interweaving constitutes, as Artin has shown, a group. 
The elements of that group are defined in the following manner. 
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Let us number at any place of the braid the lines from left to 
right, 1, 2,.... If the first line crosses the second over it, as in 
the upper crossing of the Figure 8a, we denote this process by 0, 
If the first line crosses the second underneath, as in the Mseca 
crossing in Figure 8a, we denote this process by oF ", Clearly 
a7} is the inverse of o,, for, as seen from Figure 8a, % successive 
application of the two, aye results in the reestablishment of a 
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Situation which is homotop to the initial situation. It is seen from 
Figure 8a that a mere homotopic deformation reduces 8a to just a 
set of parallel lines which do not cross at all. In a similar manner 
the relation between the second and the third line are denoted by 
o, and Os etc. If there are n lines, the number of o’s is n—1, 
the last being o,, which describes the relation between the 
(n—1)st and nth lines. The o’s constitute the generators of a 
group, the braid group, and every braid can be described as a 
power of proper products of o’s. Thus what Artin refers to gra- 
ciously as the ordinary ladies’ braid is given by 


CA o hase (7) 
That this is so is intuitively clear from Figure 8b. 


Figure 8c represents a braid; the group-theoretical expression 
for which is 


(9)0,0,2)", (8) 
while Figure 8d represents a braid 
(o,0,0,)"". (9) 


The two examples show that the braid group is noncommutative. It 
is readily seen that braids of any complexity can thus be repre- 
sented. 

We shall considerably simplify the examples to follow by con- 
sidering a special case, discussed by Artin, namely, when it does 
not matter whether one line is above or below the other. Then the 
0,8 are the same as its. The braid of Figure 8b now becomes 

(0,05)” (10) 
and is represented in Figure 9. 

Artin considers all the lines of a braid as indistinguishable, 
which is the case when they are just abstract lines. If, however, 
each line represents the world line of a particle, then if the par- 
ticles are distinguishable physically, the lines must be distin- 
guishable also. In a drawing we may represent them in different 
colors, or by full, broken, and dotted lines, etc. 

A braid may be deformed homotopically and the metric relations 
in it, such as the actual distance along a line between points of 
intersection, will change. But the topological properties expressed 
by (7) —(10) remain invariant. It is, therefore, rather natural to 
consider the relational properties of the organism as corresponding 
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to the topological properties of the four-dimensional braid of world 
lines of the particles, of which the organism is composed. A braid 
being a special kind of a knot, it may well turn out to be literally 
true that life is a knotty problem! 

It is, however, readily seen that the concept of an organism as 
a simple braid of world lines is utterly inadequate. An organism 
undergoes continuous catabolism and anabolism. While the chemi- 
cal constitution of an organism may remain relatively constant, no 
constituent atom or molecule remains in the organism for any 
length of time. Some molecules are continuously lost and are just 
as continuously replaced by other similar molecules. In terms of 
a braid, we have a braid which is constantly unwoven and rewoven 
from other threads, something like a ‘‘Belgian lace.’’ This, at 
first glance, seems to complicate the topological problem tremen- 
dously. Actually it does not. 

At any point of a braid we may ‘‘pull out’’ a thread or line homo- 
topically without disturbing the topological properties of the braid. 
Thus the braids of Figures 8b and 10a are topologically the same, 
being homeomorph. We may make the loop larger and larger, and 
in the limit move the point a (Figure 10a) into infinity, still pre- 
serving the topological properties. In the limit we obtain the situ- 
ation shown in Figure 10b, where the braid loses a line (/,) and 
gains one (/,). This can be done to any line of the braid at any 
point, except a point of intersection. In a finite four-dimensional 
world even the removal of the point a to infinity is not necessary. 
A mere removal to a sufficient distance is enough. 

Thus in our following discussion we may disregard the constant 
unweaving and reweaving of the braid and consider the classical 
form of Artin. 

An organism usually grows when it assimilates more molecules 
than it loses to the environment. In terms of a braid this means 
that more lines are woven in than are unwoven. In the situation 
shown in Figure 10, for each line lost there is one gained. We 
need, therefore, consider only the weaving in of the excess lines. 
This can be done in many different ways. Since, as we said above, 
the lines now are all physically distinguishable, we may, for ex- 
ample, postulate that every line which is at any given position on 
the outside of the braid ‘‘catches’’ a similar line unless there are 
already two lines of that kind in the braid. Such a ‘‘éaught’’ line 
then remains adjacent to the other, making a ‘‘double thread’’ of 
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identical lines until all lines are duplicated. The process is il- 
lustrated for three lines in Figure 11. At a the full and the dotted 
lines each ‘‘catch’’ a similar one; at b the broken line which here 
becomes an outside line, ‘‘catches’’ a similar line. After those 
two steps, all lines are duplicated. If there are more than three 
lines, there will be more than two steps, and the number of steps, 
in general, depends on the structure of the braid, that is, on the 
structure of its group. Each time a new line is woven in, in this 


: 


FIGURE 9. 


/ 


FIGURE 10b. FIGURE 11. 


way, a new generating element o is introduced into the braid group. 
When each line is duplicated, we must make some further assump- 
tion as to what happens. For example, we may assume that, from 
that moment on, the braid is determined by the group of all 6! 
permutations of the six threads of which it now consists, the 
permutations being obtained by a prescribed sequence of trans- 
positions. After a definite number of steps the situation will be 
reached in which we shall have the sequence: full, broken, dotted; 


54 N. RASHEVSKY 


full, broken, dotted. In other words, the original braid 1s now du- 
plicated in two similar ones. From this moment the interweaving 
may continue only amongst the first three and the second three 
lines, and thus the process begins all over again with two identi- 
cal braids instead of one. We have here a topological model of 
the duplication and continued multiplication of an organism through 
assimilation. 

Instead of the above assumptions about the method of ‘‘catching’”’ 
and weaving in of new lines, we could have made an infinite num- 
ber of different ones which lead to an eventual duplication of the 
braid. The important part is that each set of assumptions leads, 
if translated into biological terms, to different verifiable conclu- 
sions. The number of steps between the beginning of the process 
and the duplication of every line depends, as we have said, on the 
assumption made. So does the number of steps between stage c of 
Figure 11, and the stage at which the actual multiplication begins. 
But each step would correspond to some biological process, since 
each step represents an intersection of two world lines, a space 
time coincidence, that is a physically observable phenomenon. 
Thus the number of physically observable processes between the 
different stages of the life cycle of an organism are predicted by 
the topological assumptions made. 

Instead of assuming that the interweaving in separate groups of 
three begins only when the situation 


full, broken, dotted; full, broken, dotted (11) 


is reached, we may consider the case in which it begins at the 
stage 


full, full, dotted; broken, broken, dotted. (12) 


This stage will also eventually be reached after a definite num- 
ber of steps. We have now a case which biologically corresponds 
to an unequal division with differentiation. One of the new braids 
has now an excess of full lines, the other an excess of broken 
lines. We may now consider braids of second order, in which the 
role of lines is played by the braids of the first order. We may, in 
particular, consider the case in which a braid with an excess of 
full lines interweaves with one that has an excess of broken lines, 
in the same way as the full line interweaves with the broken line 
in a first order braid. In the second-order braid, which may repre- 
sent a multicellular organism, some of the individual first-order 
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braids will stand to each other in the same relation as their cor- 
responding components. (lines) stand to each other in the first- 
order braid. But this is just the kind of relation which is found in 
biology (Rashevsky, 1955c). 

In a braid of the second order there may also be a sort of divi- 
sion of functions. In two dimensions only two first-order braids 
will be ‘‘on the outside.’’ Only those two can ‘‘catch’’ new lines 
from outside, and, as in the case of Figure 11, those new lines may 
‘‘move inward’’ by a process of interweaving. When such a line 
has moved completely inward so as to become adjacent to an in- 
ward second-order braid, only then can the line be woven into that 
inner second-order braid. 

Thus the outer second-order braids may be said, to use biologi- 
cal terminology, to specialize in ingesting outside lines and to 
pass them along to the inner braids. 

Again we repeat that all the above examples are not to be taken 
seriously in any way. They are used only to illustrate a point 
which is undoubtedly very suggestive for different possible direc- 
tions of future research both in mathematical biology and in pure 
mathematics. The A-—spaces seem to be new in topology; the 
theory of braids is an established branch of the latter. As re- 
marked on p. 40, we cannot be sure that some already well de- 
veloped branch of topology does not carry in it the solution of the 
problem of geometrization of biology. 

Should such a geometrization, as anticipated in this paper, ever 
become a reality, it will be the greatest triumph of geometry and 
of pure mathematics in general. The doors to every science can 
then carry the inscription which appeared over Plato’s Academy: 


MndetG oyewNeTelnoG ElorTo. 


This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of the University of Chicago. 
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Equations are obtained which relate three important characteristics of 
a homogeneous biological aerosol; (1) the over-all viable survival; (2) 
the distribution of ages (times of suspension) of surviving cells; and (3) 
the dependence of the death rate of individual cells upon their age. The 
general solution and specific solutions are discussed and allow the de- 
termination of the kinetics of death of microorganisms to be determined 
from the experimentally measured viable cell concentration in the aerosol 
as a function of time. 


I. Introduction. The usual procedure in studying the behavior of 
a bacterial cloud is to spray a liquid suspension or dust a dry pow- 
der into a closed chamber for a finite period in order to build up 
the desired concentration; then stop generating the aerosol and 
measure the concentration of viable organisms in the air as a func- 
tion of time. Since the cloud contains organisms that have been 
suspended in the air for various lengths of time, defined as age, 
the biological decay of the whole cloud will not in general corre- 
spond to that of a cloud of organisms of the same age, unless the 
probability of death is independent of age. In order to understand 
the effect of the aerosol environment of the bacterial cell it is 
necessary to deduce the dependence of the death rate on age from 
the experimentally measured viable cell concentrations. It is the 
purpose of this paper to describe a method for doing this. 

Since the water content of droplets or particles of an aerosol 
changes drastically during its equilibration with the gas phase 
one believes intuitively that the probability of death must depend 
on age. This belief is supported by the typical nonexponential 
viable decay curves reported by Dunklin and Puck (1948) and Ferry 
and Maple (1954). 

In this report equations will be obtained which relate three im- 
portant characteristics of a homogeneous biological cloud: (1) the 
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over-all viable survival; (2) the distribution of ages of surviving 
cells; and (3) the dependence of the death rate of individual cells 
upon their age. The equations can be applied in cases where: (1) 
the aerosol generation rate varies with time; and (2) there exists 
initially any arbitrary but known age distribution of viable cells, 

Il. The fundamental equations. The number of viable cells of 
any particular age in a perfectly mixed aerosol of uniform particles 
will change because of aging, mechanical loss, and death. If 
f(t, g) is defined as the number of viable cells at the time ¢@, per 
unit of age g, f can be represented as a surface in the right handed 
t,g,f cartesian coordinate system. 

The number of viable cells f(t,g)dg in a small age range dg, at 
any time ¢, will decrease by two processes: (1) physical elimina- 
tion and death due to all causes and (2) cells will become older 
and move into the next age group. In the first process the instan- 
taneous probability of loss per unit of time, K(g), is assumed to be 
a function of age only, while for the second process the net rate of 
decrease due to aging alone is seen to be Afdg)/dg. The aerosol 
generator can only add cells of zero age and will be involved in 
the boundary conditions. Older cells can be added by the mixing 
of two aerosols, and this situation will be discussed later. 

The net rate of increase in the number of viable cells between 
the ages g and g + dg at time ¢ will be the negative sum of the 
rates of decrease due to aging and to loss: 


Hfdg) _ fdg) 


iy bale aee K(g)fdg (1) 


or 
Cie el 
i (2) 
This partial differential equation will be taken as an accurate 
description of the cloud. Its solution under the necessary boundary 
conditions will yield a surface, f(t,g) representing the age distribu- 
tion at any time. To obtain the number of viable cells suspended 
in the aerosol as a function of time it.is necessary to integrate 
f(t,9) over all possible ages. 
A general solution of equation (2) can be obtained by the method 
of Lagrange which yields: 


flt,9) = Y(t - g)eSKlelde , (3) 


The arbitrary function, w of (t-g) must be chosen so as to satisfy 
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the boundary conditions required by the physical problem. If the 
surface, f(¢,g), is known to pass through a specified curve, W, in 
general, will be defined. 

As pointed out above, the rate of loss of cells from an age group 
is simply /, so that if R(¢) is the rate at which the aerosol is being 
generated it follows from equation (3) that 


R(t) = f (2,0) = Ye, (4) 
where (0) is the value at g = 0 of 
ota) = / Ka)dg (5 


FIGURE 1 


This defines wW so that the age distribution is given by 


f(t,g) = R(t—g)e BOs (6) 


where R(t—g) is the same function of (¢—g) that F(t) was of ¢. 

As a simple example consider the age distribution for F and Kk 
constant as is often dssumed. The age distribution is plotted in 
Figure 1 for R(t) = R when 0 <¢< ¢, and zero elsewhere. 

If instead of a completely specified /(¢) an initial age distribu- 
tion, G(g) at ¢ = Ois given, the age distribution at a later time will 
be 


4 
(tg) = G(g—te Je K(g)de (1) 
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The regions of applicability of equations (6) and (7) are illustrated 
in Figure 2 in which it is assumed that & = 0 for ¢ < 0 and G =00 
for g < 0. Almost all practical applications can be made with equa- 
tion (6) and most of the following discussion will be confined to 
such cases. 

The number M#) of viable cells suspended in the aerosol at any 
time will be given by the integral of f(¢,g) over all possible ages, 
but in practice we can assume that f = 0 for 0 > g > ¢ so that 


N(t) = fr R(t— geo Kee ag, (8) 
0 


EQUATION (6) EQUATION (7) 


t=9 
FIGURE 2 


This can be considered as the basic equation describing the 
survival of microorganisms in a well mixed cloud, and it is neces- 
sary to derive K(g) from the experimentally determined M(t). _ 

Ill. General and special solutions. Let 


_é 
H(g) = ¢ Kees (9) 
so that if H(g) can be determined K(g) follows from 
H’(g) 
K(g) = -——. 10 
H(g) oe 


With this substitution, equation (8) becomes 


Ne = [Ree gH(ghdeg , (11) 
0 
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which must be solved for H(g). This is a Volterra integral equation 
of the first kind with kernel R(¢-g) and may be solved in a stand- 
ard manner using the techniques of the operational calculus. After 
taking the Laplace transform of both sides of equation (11) and ap- 
plying the Faltung theorem it is seen that 


= 
r(p) |’ 


where L~' indicates the inverse Laplace transform and n(p) and 
r(p) are the direct Laplace transforms of M(t) and R(t) respectively. 
The desired function H(g) can be obtained by substituting g for ¢ in 
H(t) calculated from (12). 

This general solution is of limited value in practice, since an 
analytical expression of the experimental survival curve, M(t), is 
usually not available even if R(¢) is carefully controlled. However, 
several special cases of particular interest will be discussed even 
though their results may sometimes be obtained from physical rea- 
soning without recourse to the general solution. 

Case I: Constant spraying rate. In this case 


H(t) = L™* (12) 


N(t) = R [unas ; (13) 
0 


where FR is a constant. From equation (12) or by direct differentia- 
tion of (13) 


N’(Z) 
H(t) = 14 
Weer (14) 
so the solution is given by equation (10) 
N**(g) 
k(g) =- —~ - (15) 
(9) Ng) 


Case II: Constant spraying turned off at t= a. In this case R(t) 
is a constant, R, for 0 < ¢< a and zero elsewhere so that R(¢~g) 
in equation (11) is zero except in the region ¢- a < g < ¢ where it 
is R. Therefore : 

t 
N(t) = R [ eee (16) 


t-a 


This reduces to the solution found in Case I for ¢< a, and for ¢>a@ 
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it is more convenient to consider equation (16) directly than to use 


the general solution. 
Assume that N(é) in the region ¢>a@ may be expanded in the 


series: 


N(t) = ) Aten ta (17) 


and that for g > 0 
H(g) = oe Bae. (18) 


where A, is positive and real. Then we have from equation (16) 
that 


t 
re Ae at = 2. RB, fi e n8dg 
B ta (19) 
= yy [eB,(e2—1)/2,] e Ant 
Therefore, 
B= AX /R(e*~1) , (20) 


so we obtain 


Tak 
H(g) = = a Be 1: (21) 
(9) R n ehnt_1 


which leads to the dependence of the loss constant on age: 


Ax, 
K(g) PCa ae seetiia (22) 
ci Aan oe \né 
2, (1) 


e hné 


The practical importance of this solution can be illustrated by 
means of data taken from Dunklin and Puck’s paper. The dashed 
line and open circles in Figure 3 are replotted from the published 
curves so that the number of viable cells is given as a function of 
the time after beginning the spray. The solid line represents N(é) 
fitted roughly by graphical methods using only two terms and as- 
Suming that spraying was constant for one minute: 


N(é) = 1100e7% 8% 4 900e79091! | for ¢ > 1 min. (23) 
Using equation (22) the death rate was calculated and plotted as 
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the solid line in Figure 4. The dashed line represents the ap- 
parent K(g) as is usually computed from the slopes of the two seg- 
ments of the decay curves. 

Case Ill: Impulsive spraying. If the cloud is generated by spray- 
ing rapidly for a very short time one believes intuitively that since 
all the cells are of about the same age the death rate of individual 
cells will be given accurately by the slope of the log-survival ver- 
sus time curve of the cloud. However, this relation can also be 
derived from the general solution by specifying a spraying rate, 
R(t), that is zero except in the region 0 < ¢< a where itis C/a. C 
is a constant and ais the short spraying time interval. In the limit 
as a—>0, R(t)-—>Cé&(t) and C—»WNa), the total number sprayed 
during the brief interval. 
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Since L4(2) = p we have that r(p) = Cp and, hence, 


HE os) = N(2)/C (24) 
Cp 
from equation (12) so that 
N’(g) 
==— e 25 
K(g) NO) (25) 


In order to determine the conditions under which an actual spray 
function can be approximated by 6-function or impulsive spraying 
consider the solution for Case II. If A,a<< 1 for all A, in the ex- 
pansion given by equation (21), H(g) is given approximately by 


¢ 1 
a. \aiiee ete Nile 26 
pepe = ne ye (26) 


Since Ra corresponds to Cin equation (24) itis seen that the spray- 
ing can be considered impulsive if it is limited to a time interval 
that is much smaller than the smallest 1/A,, where the \, are the 
‘‘decay constants’? in the expansion of M(¢). It should be noted 
that the smallest 1/A, found in the example discussed in Case II 
was 5.5 min and that the spraying time of one minute is not negli- 
gible compared to this; so that impulsive spraying could not be 
assumed. 

IV. Conclusion. Any study of the kinetics of killing of micro- 
organisms must take into account the possible effects of an age 
distribution due to the finite time required to expose the sample. 
As pointed out by the example in Case II, this is particularly im- 
portant in the study of aerosols where the environment of each cell 
changes very rapidly with time. The method proposed here, for ex- 
ample, showed an early death rate more than twice as great as that 
indicated by the method used by Dunklin and Puck. From a knowl- 
edge of the importance of the age distribution one can design ex- 
periments that readily give a valid picture of the dependence of 
the probability of death on the time of exposure. 
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A graph G may have more than one point-base B.. In a primordial 
graph P of Rashevsky’s (1954) Transformation 7, some of the point- 
bases may consist of nonspecialized points only, and some other point- 
bases: may contain specialized points. In this case, Rashevsky’s The- 
orem (1955a) on point-bases may not hold. The Theorem is certainly 
true if ali point-bases of P consist of nonspecialized points. A rigor- 
ous proof is given. Some results are derived for the more general case, 
when point-bases include both kinds of points. A general expression 
for the point-base ratio of the transformed graph P(7) is obtained. It 
is shown that with some biologically plausible assumptions Rashevsky’s 
interpretation of the point-base ratio and his conclusions are still true. 
A few simple Theorems on point-bases of graphs are included in this 
work. 


Recently N. Rashevsky (1954) has presented a transformation, 
leading from a relatively simple primordial graph P to a more com- 
plicated graph 7T(P), as a possible theoretical explanation of re- 
lations existing between biological organisms. 

Underlying Rashevsky’s work is the general principle that all 
the topological spaces or complexes by which different organisms 
are represented should be obtained from one or, at most, a few 
primordial spaces or complexes by the same transformation. For 
further comments on the choice of this transformation—clearly a 
very important point in constructing the theory—the reader is re- 
ferred to a later paper by N. Rashevsky (1955b). In this note, how- 
ever, we shall limit ourselves to the original transformation T 
given in the 1954 paper. 

It is assumed that the topological complexes describing rela- 
tions within an organism, or possibly between an organism and its 
surroundings or between members of a society, to which the trans- 
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formation T is applied, are finite, plane, directed graphs. As a 
result of the transformation we obtain again finite, plane, directed 
graphs. . 

In a subsequent paper (1955a), N. Rashevsky has examined in 
particular the problem of finding the relation between the poznt- 
bases of the primordial and the transformed graph. His result is 
embodied in Theorem 2, and the Corollary to Theorem 2 of which 
he has also given an interesting biological interpretation (/dzd., 
pp. 124, 125). It appears, however, that Rashevsky’s Theorem re- 
quires a more precise formulation because, as is made clear in 
Koénig’s book (1936), a graph P may have more than one point-base 
Bp. We shall show in this paper that the strictly correct version 
of Rashevsky’s result is the following: The point-base ratio rt, of 
the transformed graph T(P) is greater than the point-base ratio 
tp of the primordial graph P, if among all the possible point-bases 
By, of P there is none which includes specialized points. 

Before proceeding any further we recall here the definition and 
properties of a point-base in a directed graph G. By a path of the 
graph we mean a sequence of adjacent lines; a path in which all 
the lines are directed in the same manner will be called a way. In 
general for our purposes a way is sufficiently specified by indi- 
cating its end-points; a way whose end-points coincide will be de- 
noted as a cycle (or directed circuit).* If p, qs r,...are the points 
of the graph G we construct the sets II, Iq, Il,,...as follows: II, 
consists of the point p and all the points of G which can be reached 
by a way from p, and similarly for I,, Il,, etc. Performing this pro- 
cedure for all the points of the graph, we shall find a certain num- 
ber of these sets, say Il,, Il,, l,,,...ete., which are not proper 
subsets of any of the others. Such a special set, like for instance 
II,, is called a fundamental set (Grundmenge) having the points as 
its source (Quelle). A fundamental set may have more than one 
source, e.g., II, = Il,,, but indicating one of them will specify the 
fundamental set completely. On the other hand, two different fun- 
damental sets may have some points in common, but these common 
points cannot be sources. Choosing one source from each funda- 
mental set of the graph G we obtain a point-base B,, of G, which 
clearly has the following two properties: (a) any point of G can be 


* Or “‘uniformly directed cycle” in Rashevsky’s terminology (loc. cit., 
1955b, p.. 210). 
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reached by a way from a point of B,. (b) there is no way connect 
ing any two points of B,. 

Conversely, the point-base may be defined by the above two 
properties (K6énig, loc. cit., pp. 88 ff.). A graph G can, and in 
general will, therefore, have more than one point-base; all of them, 
however, have the same number of points Po, namely, as many 
points as there are fundamental sets in the graph. The quantity 
Po Will be called the point-base number. In the following we shall 
often say that a point is ‘‘a source of a graph’’ meaning that it is 
a source of a fundamental set in the graph. It is clear that each 
point of a graph belongs to at least one fundamental set (/dzd., 
p. 88., Theorem 2). 

Returning now to the primordial graph P in Rashevsky’s trans- 
formation, we shall have to consider its fundamental sets; the 
points of P are thus subdivided into sources and non-sources. But 
it will be recalled that Rashevsky also distinguishes between 
specialized and nonspecialized points of P. A priori there is no 
reason why a fundamental set should contain only nonspecialized 
points as sources. In the most general case, we can, in fact, sub- 
divide the fundamental sets into three categories: (a) those whose 
sources consist only of specialized points; (0) those whose sources 
consist only of nonspecialized points; (c) those whose sources in- 
clude both specialized and nonspecialized points. 

In P, let p’, p’’, and p’’’ be the number of fundamental sets of 
category (a), (6), and (c), respectively. Then any point-base of P 
has 

Gao i le a (1) 
points, and the point-base ratio of P is defined as 
= eed eee Te (2) 


id 
P 
X 


a, being the total number of points in the graph P. In the graph 
of Figure 1 we have, for example, p’ = p’’ = p’” = 1. 
As we have already pointed out, Rashevsky’s Theorem on point- 
bases is certainly true if 
p’ = pie = Os (3) 
Then 
on see (4) 


where r,, denotes the point-base ratio of the transformed graph. 
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FIGURE 1, First fundamental set consists of the points 1, 2, 3, 4 with 
point 1 as source (p’ = 1). Second fundamental set consists of the points 
4, 5, 6 with point 5 as source (p’” = 1). Third fundamental set includes 
the points 3, 4, 7, 8, 9 with points 7, 8, and 9 as sources (p’” = 1). 


If, however, condition (3) is not satisfied, we may have r,, > rp 
or 7,7, <7». To show this, we shall consider four very simple ex- 
amples with p’’ = 0. 

We take primordial graphs P with 4, =5 points, of which n =3 
are specialized (Figures 2a, 3, 4, 5). For the parameter m (Ra- 
shevsky, loc. cit., 1954) we choose the value 2. We recall that 
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the transformation T is determined by six steps, T,, ore ie 
Referring to Figure 2a, we shall specify T, as follows: the first 
component primordial graph (¢ = 1 in Rashevsky’s notation) spe- 
cializes in r and s (n, = 2), and loses t(n—n, =1); the second 
component primordial (7 = 2 = m) specializes in ¢ (n, = 1) and loses 
the two points r and s (n—n, = 2). As for T’,, we shall attach one 
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FIGURE 3 
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FIGURE 4 


¢ NON-SPECIALIZED POINTS 


O SPECIALIZED POINTS 
FIGURE 5 


subsidiary point to r and none to s (for <= 1); for ¢ = 2 we shall 
attach two subsidiary points to ¢. The same rules 7, and T, apply 
also to the graphs of Figures 3, 4, and 5. The residual graph con- 
sists simply of the line p — gq. 
Figure 2b shows the transformed graph obtained from that of Fig- 
ure 2a. The result of 7, (cf. Rashevsky, loc. czé. 1954) can be 
visualized as giving the primordial graph P in m (= 2) layers with 
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each specialized point fused into one; in Figure 2b the first layer 
consists of the points 7, s, ¢, p‘?, g°?; the second layer of the 
points r, s, ¢, p,q, and their connecting solid lines. After 
T,, we have added the three subsidiary points 7,, ¢,, and ¢, with 
their proper connections (dotted lines in Figure 2b), and finally by 
T,, we add the n(m —1)=3 residual graphs (1), (2), and (3) with 
their connecting lines (a), (8), and (y) (broken lines in Figure 2b). 

By these same rules the reader will easily be able to construct 
the transformed graphs of the primordials shown in Figures 3, 4, 
and 5. 

The point-base ratio of the primordial graph is equal to . in all 
our four examples; the transformed graphs have 16 points and their 
point-base ratio varies between aa and +. These results are col- 


lected in Table 1. 


TABLE 1. 


Primordial p’ of 
Graph of Primordial 


i 
5 
i 
5 
1 
5 
1 
5 


In particular we see that our second and fourth examples (Fig- 
ures 3 and 5) contradict Rashevsky’s Theorem as originally stated 
since for them we can choose a point-base B, consisting only of 
nonspecialized points, and yetr,, <r). 

It was pointed out by Rashevsky (loc. cit., 1954, 1955b) that 
cycles of graphs may play an important role in topological biology; 
for this reason it is worth-while stating two Theorems on point- 
bases and cycles; they are so simple as to require no proof, and 
yet will be found very useful in the following. 

Theorem 1: If p and q are two different sources of the same fun- 
damental set, there exists a way from p to q (including possibly 
other points of the graph) and also a way from g to 9, i.e., p and 
g are two points belonging to the same cycle. 

Theorem 2: If it is known that the point p is a source of a fun- 
damental set Il, and there exists a way from another point g of the 
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graph ¢o p, then qg is also a source of the same fundamental set 

II,, 1.e., we have Il, = II, (cf. Kénig, doc. cit., p. 89, Theorem 3). 

By Theorem 1, p and q then belong to the same cycle. All other 

points, if any, on this cycle are also sources of II, 

We recall at this point that Rashevsky’s proof of his Theorem 2 
was based on a preceding Lemma (Joc. cit. 1955a, p. 120). In the 
following let K, p, B, K’, B’, w, A, and K have the same mean- 
ing as in Rashevsky’s paper. We shall show that the Lemma as 
stated there is certainly true only if the point p which we remove 
from the graph K is not a source of K. Proof: Any point which is a 
source of Kk is also a source of kK’, This could be derived in the 
same manner as Theorem 3 below. We shall, however, give a 
slightly different proof which follows more closely, and completes, 
the one given by Rashevsky. 

1, It may be that all the points of K (except p) can be reached from 
points of B by ways that do not pass through p, in which case 
there is nothing to prove. 

2. Otherwise, the point-base of Kk’ consists of the points of B and 
those of B,. The first point-base property is obvious. As for 
the second point-base criterion we notice that 
a) there. can be no way connecting any two points of B or any 

two points of B|. 

b) having removed p, there can be no way from a point of B to 

apointofB. 

c) assume that there were a way 4, not passing through p, from 

a point r of B to a point ¢ of B. Now, all the points of 
A and of B, can be reached by ways starting at the point 
p. Therefore, no point of B, could have been a source of 
kK, since otherwise, by Theorem 2, p itself should have been 
a source. But, by the same Theorem 2, the existence of a 
way o implies that the point r is a source of K. Hence, 
there can be no way like a, and this achieves the proof of 
the Lemma. 

If the point p is a source but does not belong to a particular 
point-base B, the Lemma may or may not be true, as shown in Fig- 
ure 6. In this case, however, it is still true that the point-base 
B’ has at least as many points as the point-base B: 


PR 2 Pes (5) 
again the proof is the same as that of Theorem 3. But inequality 
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(5) is not sufficient to show that 7, > rp. Finally, if from K we re- 
move a point p which is the only source of a fundamental set I,— 
and therefore must belong to B— the point-base B’ may have one 
point less than the point-base B (cf. Figures 7a, Tb). This will 
happen if p is the only point of the fundamental set I, which does 
not belong to any other fundamental set of K. Otherwise relation 


(5) will still hold. 
5 


e< 
or 


en 


FIGURE 6, 1) Choose points r and u as point-base B of K. Point- 
base B’ of K’ consists of the same points r and u, and additional point ¢. 
2) Choose points s and u as point-base B of K. Point-base B* of K* con- 
sists of the points r, u, and ¢ Thus, there is a point s of B which is 
not included in B’, 


The residual graph Fk is obtained from the primordial P by re- 
moving all the specialized points and their connecting lines from 
P.* Let p, g, 7,...denote the nonspecialized points of P, and p, 
@ 1, ++.the corresponding points of R. To any way of RF there 
corresponds one in P, but, of course, the converse is not true since 
a way in P may pass through a specialized point. Definition: A 
primary fundamental set II- of the residual graph R is one whose 
source p corresponds to a point p which is the source of a funda- 
mental set II, in P. By Theorem 1, any other source g of Nie corre- 
sponds to a source gof Il,. 

We now prove 

Theorem 3: The residual graph has at least p’’ + p’’’ primary fun- 
damental sets (its point-base B,, therefore, at least p’’ + p’”’ 
points). 


*N. Rashevsky (loc. cit., 1955a, p. 121) includes one subsidiary point 
with the additional residual graph. We shall, however, keep the original 
definition of residual graph (Rashevsky, loc. cit. 1954, p. 333). 
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In P, let II, be a fundamental set of the second or third category, 
and the nonspecialized point p one of its sources. Then: (a) either 
p is a source of a (primary) fundamental set He in R, or (6) p is not 
a source of Rk, In this latter case, however p belongs to at least 
one fundamental set Ie of R with source g. Hence, there exists a 
connection g —> p in R, and therefore, a way g — pin P. By 
Theorem 2, g is a source of the fundamental set II, so that Vis is a 
primary fundamental set of R. 


O-— 2 — a K 
a Ss 1 
o——_~<—-e K 


FIGURE Ta. Point-base B consists of points p and s, point-base B’ of 
8 only. 


FIGURE 7b. Point-base B has two points (p and s), but point-base B’ 
has three points (s, ¢, and wv). 


Thus to each fundamental set of the second or third category in P 
there corresponds at least one primary fundamental set of &. On 
the other hand, a primary fundamental set of # cannot correspond 
to more than one fundamental set of P. We conclude that & must 
have at least p’ + p’’’ primary fundamental sets, g.e.d. In particu- 
lar, if Il, and lies have the same meaning as above, ‘it is easy to 
see that for a fundamental set Il, of the second category in P, 
there is a one-to-one correspondence between the sources of Il, 


and those of il... 
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Figure 8a shows a very simple primordial graph P formed by only 
one fundamental set of the third category. The corresponding re- 
sidual graph #, shown in Figure 8b, consists of two primary fun- 
damental sets and also of a third fundamental set. 

So far we have been dealing mainly with situations arising when 
we remove a point from a graph. We shall now prove a Theorem 
which has to do with adding some lines between two graphs. 

Let G and K be two distinct directed graphs, as shown, for ex- 
ample, in Figure 9 (solid lines). We denote by B, a point-base of 
G, and by B, one of K, The point-base numbers shall be p, and 
Px for the graphs G and K respectively. 


v 
3 * NON - SPECIALIZED POINTS 
O SPECIALIZED POINTS 
(a) (b) 


FIGURE 8 


Theorem 4: Draw any number of direct connections between ar- 
bitrary points of K and one point s of @ (Figure 9, broken lines), 
Denote the resulting graph by L and let B, be its point-base. Fur- 
thermore, in L, the points of B, shall be subdivided into two sub- 
sets: BY consisting of those p’, points pj, p%,... Dy , which 
cannot be reached by at least one way from s, and B*. including 
the remaining p’, points of B,. Thus py = pz + p’g- 

(1) If s is not a source of G, the point-base B, has 

Pr = Po + Px = Pg t+ Px - PK (6) 
points. 

(2 a) If s is a source of G, we still have Pr =Po + Px provided 
there is no way from any of the points p’ of By, to 8. 

(2 6) Otherwise, if there is at least one of the points p’ which 
has a connection ¢o s, but none from s, P,, iS given by 


Pr =Po + Px ~1=Po + Py ~ PZ 1. (7) 
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We shall give the proof in detail for case 1 above, by showing 
that a point-base B, of L consists of all points of B, and all 
points of B%,. 

a) Any point of G, including s, can be reached from some point 
of B.; any point of K from some point of By, or BY. But each 
point of B% can, in turn, be reached from s; thus, the first point- 
base criterion is satisfied by B,. 

b) s is not a source of G; hence, by Theorem 2, there is no way 
from s to any point of B.. Also, there is no way from s to a point 
of B’.. Furthermore, though there may be connections to and from 


THE POINTS @ FORM A POINT BASE IN G AND ONE IN K 


FIGURE 9 


s and B%, there is, of course, no way in K between a point of 
Bi, and one of B%- It follows from all this that in L there can be 
no connection between any two points of B,, or to and from a point 
of B, and one of BY. Finally, if pj and p’ are two different points 
of B%., a connection from p to p* could have been established only 
by passing through s, which is impossible since there is no way 
s— Di. Thus the second point-base criterion is also satisfied 
for B,. 

If s is a source of G, the point-base B, will contain one source 
s* of the fundamental set Il, of G (s and s* may coincide). By 
arguments similar to those above, it can be seen that in case 2a, 
a point-base B, consists of all points of By and all points of BY; 
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in case 2b, to obtain B, we take the points of BY, and all those of 
B,, except s*. 

Knowing B, , equations (6) or (7) follow at once. Clearly when 
equation (7) applies, the set Bi cannot be empty (px 21). Hence, 


we have quite generally: 
Po S Pp < Pg t+ Py (8) 

From the proof given above it is easy to deduce the following 
Corollary to Theorem 4: In G, assume that s is not a source and 
that it belongs to a fundamental set mio. Then any point of m6 
which is not a source of ne) cannot be a source of the completed 
graph L. 

We are now ready to prove the Theorem enunciated at the begin- 
ning of this paper. 

Denote by P, the graph obtained after, and including, ‘the step 
T, of the transformation T. As we have seen, P, contains m pri- 
mordial graphs P® arranged in m layers and ‘‘fused together’’ at 
each specialized point. Thus P, has 


AX +(m—1)(4, —n) = ma, ~n(m—1) 


points. 
The easiest way to find a point-base B, for P, is as follows. 
3 


Construct a point-base of P containing the maximum number of 
specialized points, i.e., from each fundamental set of the third 
category choose a specialized point as source. In P, take one such 
point-base from each layer. Since the specialized points are 
‘‘fused,’’ a source of the first or third category will be represented 
only once. The two point-base properties can easily be verified 
for the point-base Bp, thus constructed; altogether it has 


Pp, = p” at mp** oe ps 


points. In Bs any specialized source of the third category can be 
replaced by a nonspecialized source taken from any one, but only 
one, of the layers. 

If now P, is the graph obtained after the step 7, of the trans- 
formation T, i.e., after addition of the subsidiary points and their 
connections, we see that Bp, is a point-base also for the graph 


P,. In fact, any one of the subsidiary points can be reached by a 
way from a point of P,; on the other hand, the lines added in going 
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from P, to P; do not give rise to a way between any two points of 


ert . Therefore, 


Pps Pp =p tap +p", (9) 
and, since P, has ma, points, 
Pp 4 ada -? 
5 re Pp 
rT S——— SS 
Ps Ine, MK o a : (29) 
It is interesting to note that the point-base ratio of P, is never 
larger than that of P. We have 
p° ye Daw po 
Saeie 
tp. <p 5 a (11) 


0 


equality holding only for p’ = p’”’ = 9. 

In P, we denote the n specialized points by s,, s,,... s, (corre- 
sponding to s, eee : pf ae of the primordial); Sy shall be the 
vth subsidiary point attached to the specialized point s. The 
index v takes on the values 0, 1, 2,... li, it being understood that 
80 is the same point as s, and some of the quantities 1. may be 
zero. Since there are n(m-—1) subsidiary points we have the 


relation 
n 
ys 1. =n(m-1). (12) 


By arguments similar to those used in the proof of Theorem 3 it 
can be shown that if the specialized point ce is not a source of 
P, the corresponding point s, cannot be a source of P,. Moreover, 
in this case, none of the subsidiary points s,, (v=1, 2,...¢,) can 


be a source of P,. This follows directly from Theorem 2 since 


there is a way from s, to s,,,. 


On the other hand, ‘t s(?) is a source of the third category in P 
(and therefore s, a source of P,), each of the points s,, will be a 
source of P For, according to Theorem 1, there is a way from 

pmgta pee an aisots way from p‘? to , p being a non- 
Seecialized source of the fundamental set le). Therefore, by the 
Sk 


rules of the transformation 7, in P, there is a way from s,,, to s,, 
so that by virtue of Theorem 2 s, , must be a source of P,. 

The last step, 7,, of the transformation T, consists in adding 
one residual graph Ry to each of the subsidiary points Siy (v 21). 
Using Theorem 4 we find that the point-base number increases by a 
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certain amount @,,, (0 < < pp) as we insert each of the graphs 
R ay Thus the eee aries of the transformed graph 7T(P) is 
given by 
AOR 
Pome el ees (13) 
In many cases the quantities #, will not depend on the index », 
and may be written 2. Then equation (13) becomes 


Pp = Pp, + a Le. . (14) 


If By is a point-base of Ry and Yip the number of points of eae that 
are nae reached by a way from 8, to R. jy We shall have See a, 
Yj, %,=Y,-1, 1, depending on conditiens set forth in Theorem ver 
ae now that 


p° = aS = 0 r (3) 
Because of our assumption (3) none of the specialized points sie 
is a source of P, and thus none of the points s,, can be a source 
of P,. Also, repeated application of the Corollary to Theorem 4 
Aiace that none of the points s, can ever be a source of the graphs 
resulting in the process of adds the various R,, s. Hence case 
1 of Theorem 4 applies throughout. Moreover, ‘it is easy to see 
that equation (14) can be used instead of (138). 

Equations (2) and (14) reduce to 


ce? 


> ae (15) 
% 
Pr = mp” + 2 Lav, (16) 


By mapping the graph T(P) onto P according to the rules given 
in Rashevsky’s paper (loc. cit. 1954, p. 334), we can show that 


of Rae & 2 p”’ (17) 


for all values of i. In fact, let By be a primary source of R.. so 
that the point p™, corresponding to 2,» will be a source of P. It 
is clear that in the graph T(P) there can be no way from s Pe to DB, 


since such a way would map onto a connection sf _+ p® p~° which, 
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by Theorem 2, would imply that s® is also a source of P. Ine- 
quality (17) follows at once because there are p’’ primary-funda- 
mental sets in F.,. 

The transformed graph 7(P) has 


Ap = MX, + Um—1) (AX) ~n) (18) 
points. Using (16), (18), (17), and (12) we get 


Pr. mp’” + n(m—1)p’’ 


r= —— > ——____________ = * (19) 
i A,  MA%,_ + Mm—1) (AX, —7) _ 
Now, it is obvious that 
ir > Tp (20) 


where, we remember, 7, is given by equation (15). Indeed, since 
m > 1, inequality (20) reduces to % > M7 — n. Thus from (19) and 
(20) we finally obtain the result 


"Tr > ‘> 


which we had set out to prove. 

For the case when p’ and p’”’ do not vanish, it seems rather dif- 
ficult to obtain some general criterion which would enable us to 
decide whether 7, is greater or smaller than 7,. From (8) we de- 


t 
rive immediately the relations 


? 


p’ ae mp” a. p*” p” EY. mp”” “ por ae n(m-1)p, 


moa, + Mm—-1) (A ~n) eal mo, + n(m — 1) (% —n) 


» (21) 


where p, is the point-base number of the residual graph (of course 
Pr < %—7)- 

Subtracting the quantity 7 = pp/(%)—n) from both sides of the 
second inequality in (21) we can see that: if 


eS tp (22) 
then 
Picities (23) 
Proof: From (21) and (22) we get 
PACER 
t,— 9 = (24) 


<a ee 6 
7 tees S 
MUm—-1) (A,—n 


0 
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Because of (11), inequality (22) certainly holds if 


tT <p « (25) 
Furthermore, (25) is satisfied whenever 
Pp <Pp > (26) 


Finally, by Theorem 3, this last inequality (26) is true if p’ = 0. 
But it follows from our remarks in connection with Figure 7 that 
(26) is likely to hold even for p’ > 0, except for some rather spe- 
cial choice of the primordial graph and its specialized points. 

In particular the relation (25) will follow if the graph P is such 
that 


p< ——— (0 +9"), (27) 
(%) — 1) 
since (27) is equivalent to 
Se ge 
on 
If 
fo —fgae 


(which, as we have seen, is very pas. it follows from (11) 


and the first half of (24) that 7. - 7 < m —%, l.e., 


Tit 
By the same method used in deriving (17) we can now find a 
lower limit to the quantities a, of (13) as follows: B, and ie are 
the point-bases of the graphs Pp and Ri. It is Mite to chee 
a base B, with the maximum number of nonspecialized points, and 
such that each point of the second or third category in 5, corre- 
sponds to a (primary) point of By 


1) Let the specialized point 0 be a source of the first category 
in P. Then in P,, 8, May or may not be a source. However, it is 
clear that if pis Ae source of R., in the graph 7(P) there can be 
no way from p to s,, since this would map onto p®—+ s®), where 
pe isa ponevecialived point of P. Thus we are either in case 1 
or in case 2a of Theorem 4 so that Ty = Yip On the other hand, 
there can be no connection from s Sy to if primary point of Bi since 
the corresponding way in P would contradict the second point-base 
property of Bp. Hence, in this case, 


a 


y= yy > pte’. (28) 
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2) Now assume that s® is not a source of P. The same argu- 
ment that led to (17) shows that in this case we have again ine- 
quality (28). 

3) Finally, let s* be a source of the third category in P. Using 
as before the principle of mapping the graph 7T(P) onto P we find: 
a) There can be no way from a point 7 of Ba to 8, unless the point 
Pp corresponds to a nonspecialized source “D® of the fundamental 
set Ilcp). b) Consider the fundamental sets of the second and third 


category in P. Each of them, except possibly Il(P) gives rise to at 
J 


* NON-SPECIALIZED POINTS 


© SPECIALIZED POINTS 
FIGURE 10a 


© SPECIALIZED AND SUBSIDIARY POINTS 
FIGURE 10b 


least one primary point of he which cannot be reached by a way 
from Sips Hence, 


Vey > po” os raed Ee (29) 
c) However, it follows easily from a and b that if case 2b of Theo- 
rem 4 is realized we must have 


- 


@, +1=y,>p" +p™ : (30) 

Thus the inequality 
yp FPe t+ po’ — A (31) 
holds quite generally for this category of subsidiary points. Fig- 
ure 10 will serve to illustrate the various possibilities that may 
arise if oe is a source of the third category. Let the primordial 


graph, or ‘a partial graph of the primordial, be as shown in Figure 
10a with the three specialized points s,, s., 8, and the nonspecial- 
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ized point p as sources of the same fundamental set. Let there be 
one subsidiary point s,, for s,, two subsidiary points s,, and s,, 
for s,, and one subsidiary point s,, for s,. These subsidiary points 
are shown in Figure 10b with their connecting lines (for simplicity 
we have drawn only one layer of P,). Starting from P,, add first 
the residual graph R,,. Case 2a of Theorem 4 applies with 
ME east 4 eben = be 

Next, let us add the residual graph R,,. We are now in case 2b of 
Theorem 4, and we have y,, > p’’ + p’”’; therefore, 


oo Uap ee eer ae 


After this we add the graphs R,, and F,,. For both of them we have 
case 1 of Theorem 4, but 


a? 


Zap = Yon 2 P' +p” , 

oY > ae gatas 
If the residual graphs were added in a different order, the values 
of the @,,’S might change; however, the final result for p, would, 


of course, always be the same. 
In conclusion we find that 


p’ nm mp” ‘. pace 4 Paige + Fie = 1) + y L (on Pe pet) 
> j j 
mM, + “Km—1) (% —n) ” (32) 

where the first sum refers to those specialized points he which 
are sources of the third category in P, and the second sum to the 
remaining specialized points of P. 

The inequality (31) certainly holds for all values of j and v. In 
many cases and, in particular, if p’’’ = 0, the primordial graph may 
be such that 


> p°° +p’ for all j and v. 


On the other hand, the condition 
trae fe (33) 
can be expressed in terms of @, the mean value of the quantities 
,,. # is defined by the equation 
ee 


n(m—1) # = pf eB a (34) 


j=i vei 
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because of (31) we have 


Se oe pel, (35) 
Inserting (34) into (13) and writing out the inequality (33) we find 
pe Gor Pet) (i - x): (36) 
n XK 
Notice that 
ataletss 


lc. (37) 


Inequality (36) is completely equivalent to (33), but it has the 
advantage that the quantity % can be obtained, at least in princi- 
ple, by inspection of the primordial and the residual graph. 

Clearly, when p’’ + p’’’ — 1 is larger than the right-hand side of 
(36) we have 7, >. Another consequence of (36), not involving 
the quantity z explicitly, is the following 


Theorem 5: If a) p- = 0 (38) 
b) n> Va, (39) 

a 
ipo 2, (40) 

n 


then 1, >t. 
Proof: From (39) we obtain 


<0, (41) 


=| Ta 


n 
% 
and from (40), 


nr 
Sp eee (42) 
5 


Thus, starting from (35), using (41) and (42) we find: 


PSP n c? cad n (aad eC??? 1 n 
Se agi atic t es UP (ine ie I eles i 
0 0 0 


which is the same as (36) with p’ = 0. 
In a similar way we can derive 
Theorem 6: If, in addition to (35), the more stringent inequality 


@>p +p’ +p +1 (43) 


is satisfied, then 7, > t. 
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Proof: (36) follows easily from (43) and (37): 


- 


73 ce? 
aaa p a p 
+ ——————q—- 


n 


- 


Z>p +0 +p" + 1s p wp 0 


3 n 
> ———— (6p +p tp) t= >=}. 
n X 
These last two Theorems, 5 and 6, are sufficient to show that if 
we restrict ourselves to transformations 7 for which r, > m, this 
does not constitute a very severe limitation provided: 


ap =0, (44) 
b) a, is sufficiently large. (45) 


Consider for example the relation (39); the larger ~,, the more (39) 
is likely to be true. As for (40), we have to remember that p”” must 
always satisfy the condition 


cer 


pe <a - np "<a, —m, (46) 
It is not hard to see that the three inequalities (39), (40), and (46) 
can easily be fulfilled if a, is large enough. Again, in view of 
(35), it is clear that the inequality (43) is likely to hold for a pri- 
mordial graph with many points (more precisely, a graph with a 
large number of fundamental sets) if p’ is not too large. In particu- 
lar, (43) may be satisfied even if some of the a) 's are <p’ + p’"-+ 
are 

Finally, we shall discuss briefly the two assumptions (44) and 
(45). The second of these is a very natural one, since, as stated 
by N. Rashevsky, even the graph of a relatively simple biological 
organism is probably very complex. As for (44), it means essen- 
tially the same as Rashevsky’s assumption that specialized points 
should not be sources of the graph. In other words, a biological 
function upon which several other subordinate biological functions 
are dependent, should, in general, not become specialized. This 
again is a rather plausible hypothesis. Exceptions may arise for 
sources of the third category. But, as we have seen, such sources 
are always included in a cycle which comprises at least one non- 
specialized point. As far as the point-base properties are con- 
cerned, all the points on the cycles are equivalent; yet, biologi- 
cally, it is natural to assume that a nonspecialized point be con- 
sidered as the origin of the cycle. For instance, in Figure 3 we 
may regard the nonspecialized point p as a ‘“‘biological point- 
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base’’ or as the “‘starting point’’ of the whole graph. The line 
s—> p which closes the cycle is thought of as having been formed 
after the connections p—> g—+r—s were already present. Thus 
the biological point-bases are a subset of all possible point-bases. 
It is true that in P, each nonspecialized point is represented m 
times; nevertheless the complexity of the final graph T(P) (caused 
by addition of the residual graphs R.,) is such that in most cases 
", > Tp 

We see, therefore, from this discussion that there is no reason 
to abandon Rashevsky’s interpretation of point-bases and that, in 
general, his conclusions still hold true. It should be emphasized 
again that these conclusions were only probabilistic and may not 
hold in some particular cases. 
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BOOK REVIEW 


GEORGES ASSIE, EDITOR. French Bibliographical Digest; 2nd Series, 
6th issue, No. 14. Science; Mathematics (Part I: Pure Mathematics). 
July 1955. 128 pp. New York: The Cultural Division of the French 
Embassy. Sent on request without charge to libraries, university de- 
partments, and scientists. 


Most ordinary citizens are inclined to think of diplomats as unneces- 
Sary nuisances. It is, therefore, a pleasant surprise to find the French 
Embassy engaged in editing this Digest, which is only one of a collection 
intended primarily to make French scientific work better known in the 
United States. 

This little volume gives the titles of approximately 165 French publi- 
cations in the field of pure mathematics; some of these are symposia on 
special subjects and include several papers. It covers the period from 
1949 to 1954 and lists books, doctoral theses, and scientific papers. 

An abstract or review summarizing the content accompanies the ti- 
tle of practically each publication. Some of these abstracts—those of 
books and major papers—are quite long; all of them are competently 
written and give a clear idea of the subject matter. In fact, most ab- 
stracts are reprinted from Mathematical Reviews, the remaining, trans- 
lated into English, are taken from the Bulletin Analytique du C.N.R.S. 

All the works listed are of, at least, university level. The publications 
are ordered by subjects as follows: I. Generalities—History—Collected 
Works (page 9); II. Logic, Theory of Sets (p. 10); IIT. Algebra (p. 12); IV. 
Algebraic Geometry (p. 18); V. Number Theory (p. 25); VI. Calculus, 
Theory of Functions (‘Classical Analysis’’) (p. 34); VII. Geometry (Ele- 
mentary Geometry—Projective Geometry—Elementary Differential Geom- 
etry—General Topology) (p. 44); VIII. Homological Algebra, Algebraic 
Topology (p. 58); IX. Lie Groups and Lie Algebras, Global Differential 
Geometry, Complex Manifolds (p.'64); X. Functional Analysis (Topologi- 
cal Vector Spaces—Measures and Distributions—Operators and Rings of 
Cperators—~Harmonic Analysis—Differential, Partial Differential and 
Integral Equations—Potential Theory) (pp. 88 to 125). 

Preceding the Bibliography there is a short introduction: ‘‘Pure Mathe- 
matics in France from 1949 to 1955’? by Professor Jean Dieudonné (Uni- 
versity of Nancy, now at Northwestern University, Evanston, Illinois). It 
is a summary of the important contributions made by French mathemati- 
cians in recent years. ‘‘The most significant developments in the period 
we are considering,’’ writes Professor Dieudonné, ‘‘have occurred in 
Algebraic Topology and Functional Analysis.’’ Links and similarities, 
if any, to work in other countries are not stressed. 

The printing and other technical details are excellent. The Digest 
represents a useful task well performed. We are looking forward to read- 
ing the second part which should be published by now; it will be devoted 
to Applied Mathematics and is likely to be of even greater interest for the 
readers of this Bulletin. 

ERNESTO TRUCCO 
The University of Chicago 
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